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INTRODUCTION: 

A  new  concept  of  flat-panel  imager  (FPI)  with  avalanche  gain  and  high  resolution  (with  50  micron 
pixel  size)  is  being  investigated  for  improving  the  imaging  performance  of  digital  mammography  at 
low  dose  and  high  spatial  frequencies,  which  are  critical  for  the  detection  of  subtle  breast  abnormalities 
and  the  development  of  digital  tomosynthesis.  The  detector  employs  an  avalanche  photoconductor  - 
amorphous  selenium  (a-Se),  which  is  called  HARP,  to  detect  and  amplify  the  optical  signal  generated 
by  an  x-ray  scintillator  with  needle  structure  -  cesium  iodide  (Csl),  and  form  a  charge  image  that  is 
read  out  by  a  high-resolution  field  emitter  array  (FEA).  We  call  the  proposed  detector  SAPHIRE 
(Scintillator-Avalanche  Photoconductor  with  High  Resolution  Emitter  readout).  Our  objective  is  to 
investigate  the  feasibility  of  SAPHIRE  and  demonstrate  its  potential  advantages  over  existing  detector 
technology. 

BODY: 

The  following  tasks  were  proposed  in  the  approved  Statement  of  Work.  The  results  of  our  investigation 
in  the  first  year  are  summarized  as  follows: 

TASK  1:  Determine  the  optimal  CsI-HARP  combination  for  mammography,  Month  1-24 

a.  Develop  a  linear  system  model  for  SAPHIRE,  so  that  the  imaging  performance  of  different 
CsI-HARP  combination  can  be  predicted.  Month  1-6 

A  cascaded  linear  system  model  was  developed  for  SAPHIRE  in  order  to  determine  whether  it  is 
advantageous  to  use  high  resolution  (HR)  type  Csl,  which  has  not  been  used  in  existing  FPI  due  to  its 
low  light  output.  Since  the  avalanche  gain  of  HARP  can  compensate  for  the  low  light  output,  HR  Csl 
can  potentially  provide  better  imaging  performance  at  high  spatial  frequencies  than  the  high  light  (HL, 
optimized  for  light  output)  Csl  used  in  existing  digital  mammography  detectors.  Experimental 
measurements  of  HR  and  HL  type  Csl  samples  with  the  same  thickness  (150  micorns  for 
mammography)  showed  that  HR  Csl  has  higher.resolution  and  less  depth  dependent  image  blur,  at  the 
cost  of  lower  light  output  (60  %  of  HL).  This  resulted  in  superior  detective  quantum  efficiency  (DQE) 
at  high  spatial  frequencies  for  the  HR  Csl.  This  work  is  described  in  detail  in  a  Med.  Phys.  paper 
(manuscript  1). 

b.  Make  a  test  system  for  experimental  measurement  of  the  imaging  performance  of  CsI- 
HARP  combination  using  a  HARP  tube  with  fiber  optic  faceplate.  Month  1-12 

We  made  progress  in  the  following  two  aspects:  1.  the  first  HARP  tube  with  FOP  faceplate  was  built, 
which  allows  us  to  put  a  Csl  in  direct  contact  with  HARP  without  loss  of  light  or  spatial  resolution  for 
investigating  the  imaging  performance  of  the  scintillator  (CsI)-HARP  combination.  The  first  FOP- 
HARP  has  a  HARP  layer  thickness  of  4  pm.  Figure  1  (a)  shows  the  microscopic  view  of  the  FOP 
target  surface  after  optical  polishing  (which  cleans  the  surface  but  does  not  remove  the  surface 
variation,  i.e.  bumps).  The  “chicken-wire”  pattern  of  FOP  is  clearly  visible.  Figure  1  (b)  shows  the  first 
optical  image  made  with  the  FOP-HARP  tube  under  an  avalanche  gain  of  20.  The  result  from  the  first 
FOP-HARP  tube  is  very  encouraging  despite  the  large  number  of  defects.  We  expect  significant 
reduction  in  the  number  of  defects  after  mechanical  polishing  is  applied  to  the  FOP  for  the  next 
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prototype.  2.  Theoretical  investigation  of  the  imaging  performance  of  the  scintillator  HARP  (SHARP) 
combination  is  performed  to  determine  the  optimal  design  parameters  and  operational  conditions  for 
the  HARP  layer.  We  found  that  a  HARP  layer  with  thickness  of  4-8  pm  will  provide  adequate 
avalanche  gain  and  quantum  efficiency  for  x-ray  quantum  noise  limited  performance  down  to  a  single 
x-ray  photon.  We  also  showed  that  a  HARP  layer  thickness  variation  of  <  2  %  is  desirable  to  maintain 
a  uniform  avalanche  gain.  Direct  x-ray  interaction  in  HARP  was  shown  not  to  affect  image  quality. 
This  work  was  presented  at  the  2005  SPIE  Medical  Imaging  Conference  (manuscript  2).  A  reprint  of 
the  paper  is  attached  to  this  report.  An  expanded  version  of  this  paper  has  been  submitted  to  Med. 
Phys.  for  review,  (manuscript  3). 


Figure  1:  (a)  Micrograph  of  the  FOP  target  surface  after  optical  polishing.  The  surface  variation  is 
estimated  to  be  0.6  pm.  (b)  Optical  image  of  a  resolution  target  made  by  the  FOP-HARP  tube.  Defects 
are  most  likely  due  to  surface  variations  on  the  FOP  that  caused  local  high  field  and  breakdown  of  the 
blocking  layer.  Mechanical  polishing  or  planarization  techniques  are  under  development. 

c.  Measure  the  x-ray  imaging  performance  of  CsI-HARP  using  Csl  of  different  types  and 
thickness,  and  determine  the  optimal  CsI-HARP  combination.  Month  13-24 

TASK  2:  Investigate  the  factors  affecting  the  spatial  resolution  of  the  FEA  readout  method,  Month  1- 
36 


a.  Develop  an  electric  field  simulation  for  the  electron  optics  of  SAPHIRE  to  calculate  the 
electron  beam  trajectory  for  FEA  readout  method.  Month  1-12 

We  performed  preliminary  theoretical  investigation  of  the  factors  affecting  the  resolution  of  the  FEA 
readout  method.  Our  results  showed  that  the  dominant  source  of  image  blur  for  the  FEA  readout 
method  is  the  spread  of  the  electron  beam.  The  magnitude  of  electron  spreading  depends  on  the 
following  factors  of  the  SAPHIRE  design:  (1)  emission  angle  of  the  electrons  from  the  FEA;  (2)  design 
of  electron-optics  for  focusing  of  electron  beam;  (3)  operating  potential  of  the  FEA;  and  (4)  the 
distance  between  the  FEA  and  the  HARP  layer.  The  results  of  our  preliminary  calculation  showed  that 
incorporating  a  mesh  electrode  half  way  between  the  FEA  and  the  HARP  target  can  reduce  the  electron 
beam  spread,  and  an  effective  pixel  size  of  50  -  100  microns  is  possible.  Shown  in  Figure  2  are  the 
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calculated  electron  beam  spread  as  a  function  of:  (a)  the  distance  between  the  gate  of  the  FEA  and  the 
HARP  target;  and  (b)  the  bias  potential  on  the  mesh  electrode.  It  shows  that  to  minimize  blur  for  a 
pixel  size  of  50  pm,  an  electron  beam  spread  of  20  pm  or  less  is  desirable.  This  translate  to  a  HARP  to 
FEA  distance  of  0.4  mm  and  a  meash  electrode  bias  potential  of  350  V.  This  short  distance  with  high 
potential  on  the  mesh  will  pose  tremendous  challenge  to  the  assembly  of  the  detector.  We  are  currently 
exploring  new  electron  optics  design  that  permits  larger  distance  between  FEA  and  HARP  through 
electron  optical  simulation.  This  work  will  be  presented  at  the  upcoming  AAPM  annual  conference. 


Distance  between  gate  and  HARP  target  (mm) 


(a)  (b) 

Figure  2  Calculated  maximum  electron  beam  spread  as  a  function  of:  (a)  distance  between  the  gate  of 
the  FEA  and  the  HARP  target,  and  (b)  the  bias  potential  for  the  mesh  electrode. 

b.  Determine  the  optimal  mesh  electrode  configuration  that  results  in  the  highest  spatial 
resolution.  Month  13-18 

c.  Refine  optimal  mesh  electrode  configuration  in  connection  with  the  noise  model.  Month 
19-21 

d.  Calculate  the  spatial  resolution  due  to  FEA  readout  for  the  prototype  detector  to  verify 
experimental  measurement.  Month  25-36 

TASK  3:  Investigate  the  factors  affecting  the  electronic  noise  and  readout  speed,  Month  1-21 

a.  Investigate  the  factors  affecting  the  readout  speed  of  SAPHIRE  and  detennine  the  readout 
time  as  a  function  of  the  number  of  ITO  stripes  and  characteristics  of  charge  amplifier 
ASIC.  Month  7-12 

b.  Investigate  the  factors  affecting  the  electronic  noise  and  develop  a  noise  model  based  on  the 
properties  of  charge  amplifiers,  the  number  of  divided  ITO  stripes,  the  readout  speed  and 
the  electron  optics  design.  Month  13-18 

c.  Optimize  the  number  of  stripes  and  the  mesh  electrode  configuration  so  that  the  best  noise 
performance  and  adequate  readout  speed  can  be  achieved.  Month  19-21 

TASK  4:  Demonstrate  the  practical  feasibility  of  SAPHIRE  by  building  a  small  area  prototype 
detector,  Month  22-36 
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a.  Construct  the  small  area  prototype  detector  according  to  the  optimal  mesh  electrode 
configuration  determined  from  TASK  3.  Month  22-24 

b.  Evaluate  the  optical  imaging  performance  of  the  prototype  detector  without  Csl  at  the  input. 
Month  25-27 

c.  Evaluate  the  x-ray  imaging  performance  of  the  prototype  detector  with  Csl  coupled  to  the 
input.  Month  28-36 


KEY  RESEARCH  ACCOMPLISHMENTS: 

•  Developed  a  cascaded  linear  system  model  for  SAPHIRE. 

•  Investigated  the  imaging  performance  of  Csl  layers  made  with  different  manufacturing 
techniques,  and  determined  that  the  Csl  optimized  for  high  resolution  (HR)  can  provide  better 
DQE  at  higher  spatial  frequencies  required  for  mammography. 

•  Investigated  the  HARP  design  parameters  that  affect  the  imaging  performance  of  CsI-HARP 
combination,  and  determined  the  avalanche  gain  required  to  produce  x-ray  quantum  noise 
limited  performance  for  digital  tomosynthesis. 

•  Performed  preliminary  investigation  of  the  factors  affecting  the  image  resolution  of  SAPHIRE. 


REPORTABLE  OUTCOME: 

Manuscripts: 

1.  Wei  Zhao,  G.  Ristic  and  J.  A.  Rowlands,  “Inherent  Imaging  Performance  of  Cesium  Iodide 
Scintillators”,  Med.  Phys.  31,  2594-2605  (2004) 

2.  Wei  Zhao,  Dan  Li,  Alla  Reznik,  B.  J.  M.  Lui,  D.  C.  Hunt,  Kenkichi  Tanioka  and  J.  A.  Rowlands, 
“Indirect  flat-panel  detectors  with  avalanche  gain:  Design  and  operation  of  the  avalanche 
photoconductor”,  Proc.  SPIE  5745,  352-360  (2005) 

3.  Wei  Zhao,  Dan  Li,  Alla  Reznik,  B.  J.  M.  Lui,  D.  C.  Hunt,  J.  A.  Rowlands,  Yuji  Ohkawa  and 
Kenkichi  Tanioka,  “Indirect  flat-panel  detector  with  avalanche  gain:  Fundamental  feasibility 
investigation  for  SHARP-AMFPI  (Scintillator  HARP  Active  Matrix  Flat  Panel  Imager)”,  Med. 
Phys.  (submitted) 

Conference  presentations: 

1.  Wei  Zhao,  K.  Tanioka,  D.  C.  Hunt  and  J.  A.  Rowlands,  “A  new  flat-panel  detector  for  low  dose 
digital  x-ray  imaging”,  IEEE  Medical  Imaging  Conference,  October,  2004,  Rome,  Italy 

CONCLUSIONS 

During  the  first  year  of  our  research  support  by  Army  BCRP,  we  have  focused  our  investigation  on  the 
understanding  of  the  factors  affecting  the  imaging  performance  of  SAPHIRE.  The  imaging 
performance  of  different  Csl  layers  demonstrated  that  Csl  optimized  for  HR  can  provide  between  DQE 
at  high  spatial  frequencies  and  is  more  suitable  for  mammography.  Investigation  using  our  cascaded 
linear  system  model  of  SAPHIRE  showed  that  the  avalanche  gain  of  a  4-8  pm  thick  HARP  layer  is 
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sufficient  to  provide  x-ray  quantum  noise  limited  performance  for  digital  tomosynthesis  (using  the  HR 
Csl  layers).  Our  preliminary  calculation  of  the  maximum  electron  beam  spreading  of  the  FEA  readout 
shows  that  a  pixel  size  of  50-  100  pm  is  possible. 

REFERENCES 
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X-ray  imaging  performance  of  structured  cesium  iodide  scintillators 

Wei  Zhaoa) 

Department  of  Radiology,  State  University  of  New  York  at  Stony  Brook,  L-4  Health  Sciences  Center, 

Stony  Brook,  New  York  11 794-8460 

Goran  Ristic  and  J.  A.  Rowlands 

Imaging  Research,  Surmybrook  and  Women 's  College  Health  Sciences  Centre,  Department  of  Medical 
Biophysics,  University  of  Toronto,  2075  Bayview  Avenue,  Toronto,  Ontario  M4N  3M5,  Canada 

(Received  16  March  2004;  revised  20  May  2004;  accepted  for  publication  23  June  2004; 
published  26  August  2004) 

Columnar  structured  cesium  iodide  (Csl)  scintillators  doped  with  Thallium  (Tl)  have  been  used 
extensively  for  indirect  x-ray  imaging  detectors.  The  purpose  of  this  paper  is  to  develop  a  meth¬ 
odology  for  systematic  investigation  of  the  inherent  imaging  performance  of  Csl  as  a  function  of 
thickness  and  design  type.  The  results  will  facilitate  the  optimization  of  Csl  layer  design  for 
different  x-ray  imaging  applications,  and  allow  validation  of  physical  models  developed  for  the 
light  channeling  process  in  columnar  Csl  layers.  Csl  samples  of  different  types  and  thicknesses 
were  obtained  from  the  same  manufacturer.  They  were  optimized  either  for  light  output  (HL)  or 
image  resolution  (HR),  and  the  thickness  ranged  between  150  and  600  microns.  During  experimen¬ 
tal  measurements,  the  Csl  samples  were  placed  in  direct  contact  with  a  high  resolution  CMOS 
optical  sensor  with  a  pixel  pitch  of  48  microns.  The  modulation  transfer  function  (MTF),  noise 
power  spectrum  (NPS),  and  detective  quantum  efficiency  (DQE)  of  the  detector  with  different  Csl 
configurations  were  measured  experimentally.  The  aperture  function  of  the  CMOS  sensor  was 
determined  separately  in  order  to  estimate  the  MTF  of  Csl  alone.  We  also  measured  the  pulse  height 
distribution  of  the  light  output  from  both  the  HL  and  HR  Csl  at  different  x-ray  energies,  from  which 
the  x-ray  quantum  efficiency,  Swank  factor  and  x-ray  conversion  gain  were  determined.  Our  results 
showed  that  the  MTF  at  5  cycles/mm  for  the  HR  type  was  50%  higher  than  for  the  HL.  However, 
the  HR  layer  produces  —36%  less  light  output.  The  Swank  factor  below  K-edge  was  0.91  and  0.93 
for  the  HR  and  HL  types,  respectively,  thus  their  DQE(0)  were  essentially  identical.  The  presam¬ 
pling  MTF  decreased  as  a  function  of  thickness  L.  The  universal  MTF,  i.e.,  MTF  plotted  as  a 
function  of  the  product  of  spatial  frequency / and  Csl  thickness  L,  increased  as  a  function  of  L.  This 
indicates  that  the  light  channeling  process  in  Csl  improved  the  MTF  of  thicker  layers  more  signifi¬ 
cantly  than  for  the  thinner  ones.  ©  2004  American  Association  of  Physicists  in  Medicine. 

[DOI:  10.1118/1.1782676] 

Key  words:  indirect  detection,  flat-panel  detectors,  digital  radiography,  cesium  iodide,  Swank 
factor 


I.  INTRODUCTION 

Indirect  flat-panel  imagers  (FPI)  have  been  in  widespread 
clinical  use  for  digital  x-ray  imaging  for  several  years.  They 
use  a  scintillator  to  convert  the  incident  x-ray  energy  to  op¬ 
tical  photons,  which  are  then  converted  to  electronic  signal 
and  read  out  with  a  large  area  integrated  circuit  with  an 
amorphous  silicon  (a-Si)  photodiode  and  a  thin  film  transis¬ 
tor  (TFT)  at  each  pixel.1  Although  the  majority  of  indirect 
FPIs  use  structured  cesium  iodide  (Csl)  layers  as  the  x-ray 
converter,  powder  phosphor  screens  are  also  being  used  be¬ 
cause  of  its  lower  cost.  Unlike  the  sodium  (Na)  doped  Csl 
used  in  x-ray  image  intensifiers  (XRI1),  the  Csl  layers  used 
in  FPI  are  exclusively  doped  with  Thalium  (Tl)  because  they 
are  less  hygroscopic  and  the  emission  spectrum  (green)  is 
better  matched  to  the  spectral  sensitivity  of  the  a-Si  photo¬ 
diodes.  In  the  remainder  of  this  paper,  Csl  refers  to  Tl  doped 
layers  unless  otherwise  specified. 

Comparative  studies  between  Csl  and  powder  phosphor 
based  FPI  have  been  performed  for  both  general 


radiographic2  and  mammographic3  imaging  applications. 
The  results  demonstrated  that  Csl  provides  better  detective 
quantum  efficiency  (DQE)  than  powder  phosphors  over  all 
spatial  frequencies.  This  is  because  the  columnar  structure  of 
Csl  helps  channel  light  in  the  forward  direction.  Hence  the 
thickness  can  be  increased  while  maintaining  a  modulation 
transfer  function  (MTF)  comparable  to  phosphors.  Although 
the  imaging  performance  of  complete  detectors  incorporating 
Csl  has  been  reported,4,5  there  is  a  lack  of  systematic  study 
of  the  imaging  performance  of  Csl  alone  as  a  function  of  its 
design  parameters,  e.g.,  thickness.  Understanding  the  rela¬ 
tionship  between  physical  parameters  of  Csl  and  its  imaging 
performance  is  paramount  for  optimization  of  indirect  flat- 
panel  detectors  for  different  x-ray  imaging  applications.  The 
reason  for  pursuing  a  fundamental  understanding  of  the  im¬ 
aging  performance  of  Csl  is  that  although  Csl  has  been 
widely  used  in  commercial  flat-panel  detectors,  the  tech¬ 
niques  for  making  the  layers  and  their  relationship  with  im¬ 
aging  properties  have  largely  remained  trade  secrets.  There  is 
a  lack  of  physical  characterization  of  Csl  layers  (without 
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TFT  arrays)  and  basic  understanding  of  its  image  formation 
process,  which  is  needed  to  determine  the  optimal  Csl  thick¬ 
ness  for  a  particular  x-ray  imaging  application.  In  addition  to 
the  different  thicknesses,  commercially  available  Csl  layers 
are  of  two  different  types:  one  optimized  for  high  resolution 
(HR)  and  the  other  for  high  light  output  (HL).6  The  overall 
imaging  property  of  Csl  is  best  measured  by  its  spatial  fre¬ 
quency  dependent  DQE,  which  is  mainly  determined  by  five 
factors:  (1)  The  quantum  detection  efficiency  (QDE),  which 
depends  on  the  thickness  of  the  layer  and  the  incident  x-ray 
energy;  (2)  the  x-ray-to-photon  conversion  gain,  which  gov¬ 
erns  the  detector  performance  at  low  dose;  (3)  the  Swank 
factor,  which  describes  the  added  noise  in  the  conversion 
process  and  is  directly  proportional  to  DQE(O);7  (4)  the  over¬ 
all  spatial  resolution  expressed  as  MTF;  and  (5)  the  depth 
dependent  blur,  which  arises  because  light  generated  at  the 
x-ray  entrance  side  of  the  scintillator  experiences  more  scat¬ 
ter  before  reaching  the  exit  surface  and  more  blur  results. 
Lubberts  showed  that  depth  dependent  blur  results  in  degra¬ 
dation  of  DQE(/)  at  high  spatial  frequencies  for  powder 
phosphors.8  Compared  to  powder  phosphors,  the  columnar 
structure  in  Csl  helps  channel  light  in  the  forward  direction 
which  may  result  in  less  depth  dependent  blur.  However,  the 
light  channeling  in  Csl  is  not  perfect  as  implied  by  a  fiber 
optic-like  model.  It  has  been  shown  that  blurring  in  Csl  is 
increased  with  thickness  just  as  in  powder  phosphor  screens,6 
albeit  not  as  much.  We  have  performed  preliminary  work  on 
developing  a  model  for  the  light  scatter/channeling  process 
in  Csl  so  that  its  depth  dependent  blur  can  be  predicted.9 
Compared  to  the  isotropic  scatter  in  powder  phosphor,  the 
structured  phosphor  layer  was  modeled  to  have  nonisotropic 
scatter  with  a  preferred  direction  normal  to  the  surface  of  the 
phosphor  layer.  We  found  that  the  MTF  of  a  structured  phos¬ 
phor  was  comparable  to  a  powder  phosphor  screen  of  half 
the  thickness.  Thus  the  established  analytical  models  for 
powder  phosphors  by  Swank10  can  be  used  for  structured 
phosphors  provided  the  effective  thickness  is  modified  ap¬ 
propriately.  We  found  reasonable  agreement  between  mod¬ 
eled  and  previously  published  experimental  measurements  of 
MTF.  Others  have  developed  Monte  Carlo  techniques  to 
simulate  the  light  channeling  process  in  Csl.11  In  order  to 
validate  these  modeling  and  simulation  techniques,  system¬ 
atic  investigation  of  the  imaging  performance  of  Csl  with 
different  thicknesses  and  light  output  is  needed.  For  our 
present  investigation,  we  obtained  Csl  samples  of  four  dif¬ 
ferent  thicknesses  from  the  same  manufacturer  (Hamamatsu 
Photonics  Corporation)  for  maximum  consistency  in  material 
properties.  Only  one  of  the  thicknesses  was  available  in  both 
HR  and  HL  types,  and  the  rest  were  of  HL  type.  In  this  paper, 
we  will  present  our  experimental  measurements  of  these  Csl 
samples.  The  imaging  properties  measured  include:  (1)  The 
mass  loading  which  determines  the  QDE;  (2)  the  Swank  fac¬ 
tor  and  x-ray  conversion  gain  obtained  using  pulse  height 
spectroscopy  (PHS);  (3)  MTF;  and  (4)  DQE  as  a  function  of 
spatial  frequency.  We  also  developed  a  simplified  linear  sys¬ 
tems  model  for  Csl  which  did  not  include  the  Lubberts  ef¬ 
fect,  i.e.,  depth  dependent  blur.  By  comparing  the  modeled 
results  with  experiments,  we  obtained  qualitative  information 
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Fig.  1.  Experimental  setup  for  the  pulse  height  spectroscopy  (PHS)  mea¬ 
surements  of  Csl  samples. 


on  the  significance  of  the  Lubberts  effect  in  different  Csl 
samples.  We  believe  these  results  will  facilitate  fundamental 
understanding  of  the  light  channeling  process  in  Csl;  aid  in 
the  validation  of  different  modeling  and  simulation  tech¬ 
niques;  and  eventually  lead  to  optimization  of  Csl  param¬ 
eters  for  different  x-ray  imaging  applications. 

II.  MATERIAL  AND  METHODS 

A.  X-ray  absorption  and  Swank  factor 

1 .  X-ray  absorption  and  determination  of  mass 
loading 

The  density  of  structured  Csl  layers  is  not  well  docu¬ 
mented,  and  it  varies  between  manufacturers.  Based  on  ex¬ 
isting  documentation,  the  packing  density  of  structured  Csl  is 
~75%.12  The  nominal  thicknesses  of  the  samples  used  in  our 
experiments  were  1 50,  300,  450,  and  600  pm.  Both  HL  and 
HR  types  were  obtained  for  1 50  pm  thick  samples.  The  Csl 
layers  were  deposited  on  fiber  optic  faceplates  (FOP),  and 
are  referred  to  as  fiber  optic  scintillator  (FOS)  plates.  The 
FOP  substrates  have  a  fiber  diameter  of  6  pm,  and  they  are 
identical  for  the  HL  and  HR  type  Csl.  The  top  surface  of 
each  Csl  layer  is  coated  with  a  protective  layer  to  seal  it  from 
moisture.6  In  order  to  accurately  determine  the  x-ray  absorp¬ 
tion  of  the  samples,  we  measured  PHS  of  the  Csl  layers  for 
incident  x-ray  photons  of  different  energies,  the  results  of 
which  were  also  used  to  determine  the  Swank  factor.  Shown 
in  Fig.  1  is  our  experimental  apparatus,  which  includes  a 
photomultiplier  tube  (PMT)  to  detect  the  light  output  of  the 
FOS  plates.  The  electronic  signal  at  the  output  of  the  PMT 
(Model  R6095,  Hamamatsu  Photonics  Corporation)  was  fed 
to  a  charge  amplifier  followed  by  a  pulse-shaping  amplifier. 
The  final  output  signal  was  fed  to  a  multi-channel  analyzer 
so  that  a  PHS  could  be  obtained.  Monoenergetic  x-rays  were 
used,  which  were  generated  by  K-fluorescence  from  selected 
elements.  Each  PHS  was  integrated  to  obtain  the  total  num¬ 
ber  of  x-rays  detected  by  a  Csl  sample.  We  also  measured  the 
PHS  of  a  4  mm  Nal  (Tl)  crystal,  which  absorbs  essentially 
all  the  incident  x-ray  photons  and  provided  an  estimate  of  the 
number  of  incident  x-ray  photons.  The  ratio  between  the 
number  of  photons  obtained  using  each  Csl  sample  and  using 
the  Nal  crystal  gave  the  QDE  rj.  The  tj  values  for  each  x-ray 
incident  energy  E  was  then  fitted  to  the  expression 
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77=  1  -  exp 


0) 


where  7 i(E)/p  is  the  mass  attenuation  coefficient  of  Csl  for 
energy  E  and  d  is  the  mass  loading  or  surface  density  of  the 
Csl  layer  given  in  g/cm2.  Since  experimental  error  with  big¬ 
ger  77  values  (with  lower  E )  results  in  a  greater  error  in  the 
determination  of  d  than  for  lower  77,  a  weighting  function  of 
l/?7  was  used  in  the  curve  fitting  procedure.  The  best  fitted  d 
value  gives  the  mass  loading  of  the  Csl  sample. 


2.  Conversion  gain  and  Swank  factor 


It  has  been  shown  previously  that  the  conversion  gain  of 
HR  type  Csl  is  less  than  that  of  the  HL  type.3  This  has  a 
significant  impact  on  the  imaging  performance  of  flat-panel 
detectors  at  low  dose.  Hence  we  focused  on  the  gain  differ¬ 
ence  between  the  HR  and  HL  type  Csl  layers.  Since  the  x-ray 
conversion  gain  in  Nal  crystals  is  known,  by  comparing  the 
PHS  of  Csl  and  Nal  at  the  same  incident  x-ray  energy  E,  we 
can  obtain  the  conversion  gain  of  the  Csl  samples.  The 
Swank  factor  A  s  can  be  obtained  from  the  PHS  measure¬ 
ments  using7,13 


Ac — 


m\ 


m0m2 


(2) 


where  m0,  m , ,  and  m2  are  the  zeroth,  first,  and  second  mo¬ 
ments  of  the  PHS,  respectively. 

The  Swank  factor  describes  the  distribution  of  the  number 
of  optical  light  photons  produced  by  x-ray  photons  for  a 
given  incident  x-ray  energy  distribution  and  characterizes  the 
noise  associated  with  the  x-ray  to  light  conversion  process. 
At  each  x-ray  energy,  the  Swank  factor  can  be  separated  into 
two  multiplicative  factors:13 


^S-^OPD^AED)  0) 

where  7OPD  describes  the  optical  component  of  As  due  to 
unequal  light  losses  from  different  parts  of  the  Csl  layer 
because  of  different  propagation  lengths,  and  7AED  describes 
the  absorbed  energy  distribution  for  the  incident  x-ray  pho¬ 
tons,  which  is  mainly  due  to  K-fluorescent  escape  and  reab¬ 
sorption.  /0PD  can  usually  be  assumed  to  be  x-ray  energy 
independent,  whereas  /AED  depends  strongly  on  incident 
x-ray  energy.13  While  we  assume  70Pd  in  our  FOS  samples 
depends  mainly  on  the  properties  of  the  Csl  layer,  the  FOP 
could  also  affect  7OPD  if  the  acceptance  angle  of  the  FOP  is 
small  and  the  angular  distribution  of  light  from  Csl  is  depth 
dependent.  Since  the  FOP  used  in  our  FOS  samples  has  a 
numerical  aperture  of  l,6  we  expect  the  effect  of  FOP  to  be 
negligible.  7AED  for  a  given  x-ray  detector  material  can  be 
calculated  numerically  by  dividing  the  thickness  of  the  ma¬ 
terial  into  many  thin  layers.14  The  K-fluorescent  photons 
generated  from  each  sub-layer  are  assumed  to  be  isotropic, 
so  that  the  probability  of  reabsorption  can  be  calculated  by 
integrating  the  number  of  reabsorbed  photons  at  each 
solid  angle  and  depth.  Then  a  simulated  PHS  can  be  gener¬ 
ated  with  the  main  peak  representing  the  total  energy  absorp¬ 
tion  in  the  material  (including  K-fluorescence  reabsorption) 
and  a  second  (lower)  peak  representing  K-fluorescent  escape. 


J 


x-rays 


Fig.  2.  Diagram  showing  the  CMOS  image  sensor  and  image  acquisition 
system  used  to  measure  the  MTF,  NPS,  and  DQE  of  Csl  layers.  The  FOP 
side  the  Csl  samples  were  placed  in  direct  contact  with  the  CMOS  sensor 
(thickness  of  different  layers  not  drawn  to  scale). 


A  Poisson  distributed  noise  can  then  be  added  to  each 
peak.  Then  the  moments  ( m0 ,  m,,  and  m2)  of  the  simulated 
PHS  are  calculated  and  7AED  determined  using  Eq.  (2). 
This  method  has  been  used  in  the  theoretical  calculation  of 
7aed  for  phosphor  screens,15  Csl  in  image  intensifiers7 
and  amorphous  selenium  x-ray  photoconductors.16  After 
7aed  *s  calculated,  we  can  obtain  an  estimate  of  70Pd 
from  the  measured  Swank  factor  As.  This  information 
is  important  for  understanding  the  optical  properties  of 
different  Csl  layers. 


3.  DQE  at  zero  spatial  frequency 

It  is  difficult  to  obtain  a  precise  measurement  of  DQE(O) 
from  the  measurement  of  spatial  frequency  dependent 
DQE(/)  (Sec.  II C).  This  is  because  the  DC  component  of 
the  noise  images  used  to  compute  NPS  is  subtracted  before 
the  Fourier  transformation  is  applied.  An  estimate  of  DQE(O) 
is  usually  obtained  by  extrapolation  of  the  DQE(/)  curve 
near  zero  spatial  frequency.  Our  independent  measurement  of 
77  and  As  provided  a  more  straightforward  method  for  deter¬ 
mining  DQE(O)  using 

DQE(0)=t7^5.  (4) 

The  DQE(O)  obtained  using  this  method  is  for  monoenergetic 
x-rays.  DQE(O)  for  a  broad  x-ray  spectrum  can  then  be  esti¬ 
mated  using  the  mean  energy  of  the  entire  spectrum. 

B.  MTF 

The  MTF  of  Csl  samples  was  measured  using  a  two- 
dimensional  array  of  CMOS  optical  sensor  with  a  MOSFET 
switch  and  a  photodiode  at  each  pixel.  The  design  and  per¬ 
formance  of  the  CMOS  sensor  were  described  in  detail 
elsewhere.17  It  has  a  pixel  pitch  of  48  /x m  and  an  active  area 
of  24.5  mmX49.0  mm.  The  active  area  of  each  photodiode  is 
43  /xmX43  /xrn,  yielding  a  pixel  fill-factor  of  ~80%.  The 
FOP  side  of  each  Csl  sample  was  placed  in  direct  contact 
with  the  surface  of  the  CMOS  sensor,  as  shown  in  Fig.  2. 
Hence  the  imaging  geometry  is  similar  to  an  indirect  flat- 
panel  detector.  The  FOP  substrate  for  each  Csl  sample  is 
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3  mm  thick.  It  absorbs  99%  of  the  radiation  that  transmits 
through  the  Csl  layer,6  hence  direct  x-ray  interaction  in  the 
CMOS  sensor  is  negligible.  The  output  of  the  CMOS  sensor 
was  fed  to  a  frame  grabber  in  a  PC.  The  presampling  MTF  of 
the  CsI-CMOS  combination  was  measured  using  the  slanted 
slit  technique.18  A  10  /rm  wide  slit  made  of  1.5  mm  thick 
tungsten  sheet  was  used  in  the  experiments.  The  slit  was 
placed  at  an  angle  ~2°  from  the  gate-  or  data-lines  of  the 
CMOS  sensor.  An  x-ray  image  of  the  slit  was  taken,  and  the 
image  was  corrected  for  offset  and  gain  nonuniformity,  as 
well  as  nonresponsive  pixels.  The  maximum  pixel  value  of 
each  image  line  (perpendicular  to  the  slit)  was  determined  so 
that  those  pixels  centered  on  the  slit  could  be  located.  Their 
location  was  then  fitted  to  a  straight  line  in  order  to  deter¬ 
mine  the  position  of  the  slit  on  the  image.  Depending  on  the 
point  of  intersection  between  the  slit  and  the  pixels,  the  im¬ 
age  data  for  each  line  of  the  image  were  placed  into 
four  different  bins.  The  image  data  from  the  four  bins  were 
then  interleaved  to  produce  a  four  times  oversampled 
line  spread  function  LSF.  Then  the  data  for  the  baseline  of 
the  LSF,  i.e.,  data  less  than  1%  of  the  maximum  value  of 
the  LSF,  were  fitted  with  a  single  exponential  decay  function. 
This  reduces  the  noise  for  subsequent  Fourier  transform  of 
the  LSF,  which  produces  the  MTF.  In  order  to  determine 
the  MTF  of  the  Csl  samples  alone,  the  slanted  slit  was 
also  imaged  directly  with  the  CMOS  sensor  (without  Csl  in 
place)  using  x-rays.  This  allows  determination  of  a  correc¬ 
tion  MTF  due  to  the  focal  spot  blur,  the  width  of  the  slit 
and  the  aperture  function  of  the  CMOS  sensor.  The  correc¬ 
tion  MTF  was  then  divided  into  the  presampling  MTF  of 
the  entire  imaging  system  to  obtain  the  MTF  due  to  Csl 
alone. 

C.  DQE  measurement 

The  2D  NPS  of  the  CsI-CMOS  imaging  system  was 
measured  as  a  function  of  x-ray  exposure.  An  RQA5  x-ray 
spectrum  was  used,  which  was  generated  with  a  Tungsten 
tube  operated  at  70  kVp  (using  a  high  frequency  generator) 
and  filtered  with  21  mm  of  aluminum  (Al).  The  measured 
HVL  of  the  beam  was  7.1  mm  of  Al.  The  linearity  of  the 
detector  was  measured  within  the  exposure  range  of 
0.14-20.9  mR,  from  which  the  x-ray  sensitivity  k0  was 
determined.  At  each  x-ray  exposure  level,  100  x-ray  images 
and  dark  frames  were  acquired.  The  average  of  the  dark 
frames  formed  the  offset  image,  which  was  subtracted 
from  the  average  of  the  100  x-ray  images  to  obtain  a  gain 
correction  table.  A  defect-free  region  of  interest  (ROI) 
containing  128X128  pixels  was  chosen  for  the  NPS  calcula¬ 
tion.  Each  of  the  100  x-ray  images  was  offset  and  gain 
corrected,  and  the  mean  of  each  noise  image  subtracted.  A 
two-dimensional  (2D)  Fourier  transform  was  then  performed 
on  each  corrected  noise  image  and  the  results  were  averaged 
to  obtain  a  2D  NPS.19  Seven  lines  on  either  side  of  the  cen¬ 
tral  axes  (not  including  the  axes)  of  the  2D  NPS  were  used  to 
compute  the  one-dimensional  (ID)  NPS.19  The  frequency  of 
each  data  point  within  this  region  was  calculated  using  / 
=  \ju2  +  v2.  Although  our  method  of  gain  correction  intro- 


Oafr  Cal  Blur 


Fig.  3.  Imaging  stages  of  the  simplified  cascaded  linear  system  model  for 
the  CsI-CMOS  combination. 

duced  a  small  amount  of  correlation  between  each  noise  im¬ 
age  and  the  gain  correction  table  (~1%),  it  is  within  the 
experimental  error  of  the  determination  of  NPS  and  is  there¬ 
fore  neglected.20 

D.  DQE  model 

We  adapted  a  linear  system  model  developed  for  indirect 
flat-panel  detectors.21  Since  the  x-ray  spectrum  used  in 
our  experiments  was  70  kVp,  we  added  the  effect  of 
K-fluorescence  reabsorption.22  The  imaging  stages  involved 
in  the  model  are  shown  in  Fig.  3.  The  path  corresponding  to 
K-fluorescence  reabsorption  included  fluorescent  photons 
from  both  Cs  and  I  atoms.  We  assumed  that  if  a 
K-fluorescence  photon  interacts  with  the  material,  the  energy 
is  completely  reabsorbed  at  the  interaction  site.  Since  the 
K-Beta  fluorescence  photon  for  Cs  is  above  the  K-edge  of  I, 
3%  and  4%  of  the  photons  above  the  K-edge  of  Cs  could 
generate  a  second  K-fluorescence  depending  on  the  thickness 
of  Csl,  however,  this  was  ignored  in  our  calculations.  The 
overall  error  caused  by  the  omission  on  the  MTF  and  DQE 
should  be  negligible.  For  the  present  paper,  we  also  ignored 
the  Lubberts  effect  in  our  linear  system  model.  By  compar¬ 
ing  modeled  results  to  experimental  measurement  of  DQE, 
we  can  obtain  a  qualitative  assessment  of  the  importance  of 
Lubberts  effect  in  different  Csl  layers. 

III.  RESULTS  AND  DISCUSSION 

A.  X-ray  absorption  and  Swank  factor 

1.  X-ray  absorption  and  determination  of  mass 
loading 

Shown  in  Fig.  4  is  the  measurement  of  QDE  77  of  the  1 50 
pm  HL  type  Csl  layer  as  a  function  of  the  incident  x-ray 
energy  E.  Plotted  in  the  same  graph  is  the  calculated  QDE 
using  the  fitted  mass  loading  d  of  the  same  Csl  layer,  which 
was  determined  as  55  mg/cm2.  The  quality  of  fitting  is  ex¬ 
cellent.  Table  I  summarizes  the  fitted  mass  loading  d  for  all 
the  Csl  layers  used  in  our  experiments.  Shown  in  the  same 
table  is  the  mass  loading  determined  from  weight  measure¬ 
ments  provided  by  the  manufacturer  for  each  Csl  sample. 
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Fig.  4.  Comparison  of  experimentally  measured  quantum  detection  effi¬ 
ciency  (QDE)  as  a  function  of  x-ray  energy  for  the  150  HR  Csl  layer  and  the 
calculated  QDE  using  the  fitted  mass  loading  55  mg/cm2. 

The  mass  loading  determined  from  our  x-ray  measurement  is 
in  good  agreement  with  the  manufacturer’s  weight  measure¬ 
ment.  Compared  to  the  150  /xm  layers,  the  mass  loading 
for  the  300  /xm  layer  is  less  than  what  is  expected  from 
its  doubled  thickness.  The  most  noticeable  fact  is  that 
the  450  and  600  /xm  layers  have  very  little  difference  in  their 
mass  loading.  This  is,  however,  within  the  manufacturers 
tolerance  of  ±10%.  If  we  assume  a  constant  density  of  3.38 
g/cm3  (i.e.,  75%  packing  density),  the  thickness  L  of  differ¬ 
ent  samples  can  be  estimated  from  the  measured  mass  load¬ 
ing  and  the  results  are  shown  in  Table  I. 

2.  Conversion  gain  and  Swank  factor 

Shown  in  Fig.  5  are  the  PHS  measured  at  £  =  31  keV 
for  the  150  /xm  HR  and  HL  Csl  layers,  as  well  as  the 
Nal  crystal.  The  relative  position  of  the  peak  of  each  PHS 
determines  the  conversion  gain  difference  between  these 
three  samples.  Since  E  is  lower  than  the  K-edge  of  Csl, 
there  is  only  one  peak  in  the  PHS  of  Csl.  The  conversion 


gain  of  the  HR  type  is  —64%  of  the  HL  type.  This  ratio 
is  higher  than  that  (—50%)  observed  in  a  previous 
study,3  which  could  simply  be  due  to  variation  between 
samples.  The  position  of  the  PHS  peak  for  the  HL  type  Csl  is 
—47%  of  that  for  the  Nal  crystal  but  before  direct  compari¬ 
son  this  has  to  be  corrected  for  the  different  spectral  outputs 
of  the  Nal  and  Csl.  The  accepted  light  yield  in  Nal 
for  gamma  rays  is  —40  photons/keV,23  which  means  that 
the  energy  required  to  generate  a  light  photon  is  W 
=  25  eV.  Shown  in  Fig.  6  are  previously  published  emission 
spectra  for  both  Nal  (Tl)  and  Csl  (Tl).24  The  peak  of  the 
spectra  is  415  nm  and  530  nm  for  Nal  and  Csl, 
respectively.23  The  quantum  efficiency  (QE)  of  the  PMT 
used  in  our  experiments25  is  also  plotted  as  a  function  of 
wavelength  in  Fig.  6.  It  shows  that  the  optical  QE  of  PMT  is 
better  matched  to  the  emission  spectrum  of  Nal.  We  calcu¬ 
lated  the  average  QE  of  the  PMT  for  each  emission  spectrum 
by  integrating  the  product  of  light  intensity  and  QE  over  the 
whole  spectrum  and  normalized  it  to  the  total  number  of 
integrated  light  photons.  Our  results  showed  that  the  average 
QE  of  the  PMT  for  the  emission  spectra  of  Nal  and  Csl  are 
20.3%  and  10.8%,  respectively.  This  information  combined 
with  the  positions  of  the  PHS  peak  for  Nal  and  Csl  in  Fig.  5 
leads  to  the  conclusion  that  the  effective  W  for  the  HL  Csl 
FOP  is  —28  eV.  Since  the  FOP  substrate  of  the  Csl  transmits 
only  —60%  of  the  light  emitted  from  Csl,6  the  W  associated 
with  the  inherent  conversion  gain  of  HL  Csl  is  —17  eV.  This 
is  in  good  agreement  with  previously  published  values  for 
columnar  Csl.26 

Shown  in  Fig.  7  are  the  PHS  of  the  150  /xm  HR  Csl 
layer  at  different  x-ray  incident  energies.  For  E  above  the 
K-edge  of  Csl,  there  are  two  peaks  on  the  PHS,  with  the 
lower  one  corresponding  to  the  x-ray  interaction  events 
with  the  escape  of  K-fluorescence.  The  moments  of  the 
PHS  were  calculated  so  that  the  Swank  factor  could  be 
determined  using  Eq.  (2).  Shown  in  Fig.  8(a)  is  the  calcu¬ 
lated  Swank  factor  As  plotted  as  a  function  of  E  for  the 
150  /xm  HR  layer  and  also  the  calculated  /AED  using  its 
measured  mass  loading  of  55  mg/cm2.  The  /AED  is  unity 
below  the  K-edge  because  there  is  no  K-fluorescence  escape. 
Above  the  K-edge,  the  values  for  /AED  increase  as  a  function 


Table  I.  Nominal  thickness,  measured  weight  and  type  of  Csl  samples  used  in  our  experiments  as  specified  by 
the  manufacturer.  The  measured  mass  loading  was  determined  from  x-ray  quantum  detection  efficiency  of  each 
sample  using  x-ray  pulse  height  spectroscopy.  The  estimated  thickness  for  each  sample  was  calculated  by 
assuming  a  constant  density  of  3.38  g/cm3,  which  means  a  packing  density  of  75%. 


Nominal  thickness 

W 

150 

150 

300 

450 

600 

Mass  loading  from 
weight 
(mg/cm2) 

50 

51 

88 

142 

181 

Type 

HR 

HL 

HL 

HL 

HL 

Mass  loading  from 
attenuation 
(mg/cm2) 

55 

57 

83 

143 

166 

Estimated  thickness 
(/zm) 

163 

169 

246 

423 

491 

^OPD 

0.88 

0.90 

0.91 

0.93 

0.95 
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Fig.  5 .  PHS  measurements  at  3 1  keV  for  the  1 50  /tm  HR  and  HL  Csi  layers, 
as  well  as  the  Nal  crystal.  ' 


Fig.  6.  The  emission  spectra  of  both  Nal  and  Csi  (Tl)  scintillators  plotted  as 
a  function  of  the  wavelength  of  light.  Plotted  in  the  same  graph  is  the 
quantum  efficiency  of  the  PMT  (Hamamatsu  Model  R6095)  (scale  shown 
with  they-axis  on  the  right)  as  a  function  of  wavelength. 


of  x-ray  energy  E.  This  is  because  the  relative  position  of  the 
escape  peak  is  closer  to  the  main  peak  as  E  increases. 
The  calculated  /AED  was  then  multiplied  by  a  constant  /0PD 
until  the  result  matched  well  with  the  measured  As-  This 
provided  a  theoretical  estimate  of  As  and  is  shown  in 
Fig.  8(a),  and  corresponds  to  a  constant  /0PD  of  0.88. 
Although  the  calculated  As  matches  the  measurements 
very  closely  above  the  K-edge,  it  underestimates  the  As 
below  the  K-edge  (^=0.91)  by  ~3%.  Considering  that 
/0PD  is  energy  independent  is  an  approximation,  this  is  a 
reasonable  agreement.  Figure  8(a)  shows  that  below  the 
K-edge,  optical  factor  /0PD  is  the  only  mechanism  for 
degradation  in  As.  K-fluorescence  escape,  i.e.,  7AED  be¬ 
comes  the  dominant  factor  for  A  s  just  above  the  K-edge,  and 
then  becomes  less  important  as  E  increases.  In  Fig.  8(b)  is 
the  measured  As  for  all  Csi  layers  used  in  our  experiments, 
calculated  from  the  PHS  measurements  using  Eq.  (2).  Plotted 
in  the  same  graph  is  a  theoretical  estimate  of  A  s  obtained  by 
multiplying  the  calculated  /AED  for  each  layer  by  a  constant. 
This  produced  an  estimate  for  70PD  for  each  Csi  layer,  and 
the  values  are  summarized  in  Table  I.  In  addition  to  the 
energy  dependence,  Fig.  8(b)  also  shows  an  increase  in  As 
for  thicker  Csi  layers.  This  is  because  of  the  increased 
probability  of  K-fluorescence  reabsorption,  which  results  in  a 
higher  IAED .  Also,  as  shown  in  Table  I,  the  estimated  fOPD 
values  increase  as  a  function  of  thickness,  and  the  70Pd  for 
all  the  HL  layers  are  higher  than  that  for  the  HR.  This  sug¬ 
gests  a  difference  in  optical  properties  between  different  Csi 
layers,  e.g.,  the  number  of  light  photons  escaping  from  dif¬ 
ferent  depths.  It  is  important  to  note  that  the  Swank  factors 
for  the  150  jum  HR  and  HL  Csi  layers  are  very  close  (0.91 
and  0.93,  respectively,  below  the  K-edge).  This  differs  from 
powder  phosphors,  where  high  resolution  is  usually  achieved 
by  replacing  the  reflective  backing  with  an  absorber,  result¬ 
ing  in  additional  gain  fluctuation  noise  and  a  much  lower 
Swank  factor  which  can  approach  0.5. 15  The  high  Swank 
factor  of  HR  Csi  suggests  that  even  if  an  absorptive  backing 
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is  used,  the  light  output  from  deep  within  the  layer  is  not 
reduced,  i.e.,  there  is  very  little  variation  in  the  amount  of 
light  emitted  per  absorbed  x-ray  photon  from  different  depth 
of  the  Csi  layer.  This  needs  to  be  considered  when  develop¬ 
ing  physical  models  for  different  types  of  Csi  layers.  Another 
possibility  is  that  both  the  HR  and  HL  Csi  layers  have  the 
same  reflective  backing,  but  the  HR  layers  have  light  absorb¬ 
ing  materials  between  the  Csi  columns  to  minimize  lateral 
spreading  of  light. 

Although  HR  Csi  has  been  developed  for  some  time,  it 
has  not  been  used  in  commercial  flat-panel  detectors  (includ¬ 
ing  mammography).  This  is  perhaps  due  to  its  lower  conver¬ 
sion  gain,  which  makes  the  FPI  more  susceptible  to  elec¬ 
tronic  noise  and  reduces  DQE  at  low  exposure  levels.3,27 
Hence  the  resolution  advantage  of  the  HR  type  Csi  has  not 
yet  been  realized  in  the  clinical  products.  Development  of 
new  indirect  flat-panel  detectors  with  avalanche  gain  can  po¬ 
tentially  compensate  for  the  gain  loss  in  HR  Csi,28  so  that  it 
can  be  used  for  high  resolution  x-ray  imaging  applications, 
e.g.,  mammography,  in  order  to  boost  the  DQE  at  high  spa¬ 
tial  frequencies. 

3.  DQE  at  zero  spatial  frequency 

From  the  experimentally  measured  rj  and  As,  DQE(0) 
for  different  Csi  samples  was  calculated  using  Eq.  (4).  In 
order  to  obtain  a  smooth  curve,  the  fitted  mass  loading 
d  shown  in  Table  I  and  the  corresponding  r)  for  different 
incident  x-ray  energy  E  were  used.  Shown  in  Fig.  9  is 
the  measured  DQE(0)  as  a  function  of  incident  x-ray  energy 
for  all  the  Csi  samples  used  in  our  experiments.  DQE(O) 
for  the  1 50  pm  HR  and  HL  layers  are  essentially  identical, 
which  is  to  be  expected  from  their  similar  thickness 
and  Swank  factor.  DQE(0)  increases  as  a  function  of  the 
thickness  due  to  increases  in  both  77  and  A  s .  Below  the 
K-edge  for  each  Csi  layer,  DQE(0)  decreases  as  a  function  of 
x-ray  energy  E  due  to  the  drop  in  77.  DQE(0)  exhibits  a 
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Fig.  7.  PHS  measured  at  nine  different  x-ray  energies  ranging  from  25  to  74  keV  for  the  150  fim  HR  Csl  layer. 


significant  jump  just  above  the  K-edge  due  to  the  increase  in 
r/,  despite  the  decrease  in  As.  However,  DQE(0)  decreases 
again  as  a  function  of  E  for  energies  farther  above  the 
K-edge. 

B.  MTF 

Shown  in  Fig.  10(a)  is  the  measured  MTF  of  different  Csl 
samples.  This  is  after  the  correction  MTF  (i.e.,  pixel  aperture 
function,  slit  width  and  the  effect  of  focal  spot  size),  which  is 
also  shown  in  Fig.  10(a)  for  reference,  has  been  removed 
from  the  presampling  MTF  of  the  CsI-CMOS  combination. 
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The  150  /rm  HR  Csl  layer  has  the  highest  MTF,  as  expected. 
The  MTF  at  5  cycles/mm  for  the  150  /mi  HR  layer  is  50% 
higher  than  that  of  the  150  fi m  HL.  The  MTF  of  the  HL  type 
Csl  layers  decreases  as  a  function  of  thickness. 

There  is  another  factor  affecting  the  thickness  dependence 
of  MTF  in  Fig.  10(a),  the  blur  caused  by  K-fluorescence 
reabsorption.  It  has  been  shown  from  the  consideration  of  A  s 
that  K-fluorescence  reabsorption  depends  on  Csl  thickness. 
Since  the  spectrum  used  in  our  experiment  (RQA5)  has  the 
majority  of  its  x-ray  energy  above  the  K-edge  of  Cs  and  I, 
K-fluorescence  plays  a  role  in  the  MTF.  In  order  to  obtain  a 
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20  30  40  60  60  70  80 

(b)  X-ray  energy  (keV) 

Fig.  8.  (a)  Swank  factor  for  the  150  /am  HR  layer:  The  solid  squares  show 
the  experimentally  measured  A s  from  the  PHS,  the  solid  line  is  the  calcu¬ 
lated  /AED ,  and  the  dashed  line  is  the  /AED  scaled  by  a  constant  of  0.88.  (b) 
Swank  factors  for  all  the  Csl  samples:  The  solid  symbols  represent  experi¬ 
mentally  determined  A  s ,  and  the  lines  represent  calculated  A  s  values  by 
scaling  the  /AED  due  to  K-fluorescence  reabsorption  by  a  constant  represent¬ 
ing  the  estimated  /opd  •  The  thickness  values  shown  in  the  label  are  the 
nominal  thickness  given  in  /xrn  for  each  layer.  The  type  of  the  Csl  layer  is 
HL  unless  specified  as  HR. 


reasonable  estimate  of  the  MTF  due  to  optical  factors  alone, 
we  computed  the  MTF  due  to  K-fluorescence  reabsorption 
for  each  Csl  layer  using  a  numerical  method  developed 
previously.29  The  resulting  MTFK  are  shown  in  Fig.  10(b). 
The  MTF  shown  in  Fig.  10(a)  were  then  divided  by  the 
MTFk  so  that  the  MTF0,  i.e.,  MTF  due  to  optical  factors 
alone  could  be  estimated.  The  results  for  MTF0  are  shown 
in  Fig.  10(c)  for  each  Csl  layer.  Since  MTFK  decreases 
in  thicker  layers  due  to  higher  probability  of  K-fluorescence 
reabsorption,  the  correction  is  more  important  in  the 
MTF0  for  thicker  layers  than  the  thinner  ones.  As  shown 
in  Fig.  10(c),  MTF0  for  the  600  /am  HL  sample  is  essentially 
the  same  as  that  for  the  450  /am  sample.  This  is  partially 
because  there  is  very  little  difference  between  their  estimated 
thicknesses.  Similarly,  the  MTF  for  the  300  /am  HL  layer 
is  very  close  to  the  150  /am  HL.  After  communicating 


Fig.  9.  Calculated  DQE(0)  from  the  measured  Swank  factor  As  and 
x-ray  quantum  detection  efficiency  p  (used  fitted  mass  loading  shown  in 
Table  I). 


the  results  to  the  manufacturer,  we  were  informed  that 
the  300  and  600  /am  HL  samples  were  manufactured  in  a 
different  deposition  system  and  had  a  slightly  thinner  colum¬ 
nar  diameter  than  the  rest  of  the  samples.  This  also  explains 
why  the  MTF  of  these  two  layers  are  so  close  to  the  next 
thinner  layer. 

As  shown  in  Fig.  10(c),  although  the  MTF  for  Csl 
layers  of  different  thickness  (HL  type)  differ  significantly  at 
medium  and  high  spatial  frequencies,  their  MTF  below  1 
cycles/mm  are  surprisingly  similar.  Our  hypothesis  is  that 
a  fraction  of  the  Csl  layer  lacks  the  columnar  structure. 
The  most  likely  place  for  this  layer  to  occur  is  at  the  begin¬ 
ning  of  deposition  (the  layer  closest  to  the  FOP  in 
our  samples).  This  is  because  the  FOP  substrate  of  our 
samples  was  not  patterned  prior  to  Csl  deposition,  and  the 
columnar  crystal  structure  only  starts  to  form  after  a  nominal 
thickness  of  a  structureless  layer  is  deposited.30  We  speculate 
this  structureless  Csl  scintillator  layer  may  act  as  a  transpar¬ 
ent  layer  (or  at  best  as  a  powder  phosphor),  both  of  which 
are  known  to  have  poorer  MTFs  than  structured  scintillators. 
This  explains  why  the  MTF  of  all  the  Csl  layers,  regardless 
of  their  thickness,  drops  rapidly  at  low  spatial  frequency. 
Given  that  each  Csl  layer  is  composed  of  a  fixed  thickness 
of  structureless  layer  in  addition  to  the  columnar  structure, 
the  fraction  of  structured  Csl  increases  as  a  function  of 
thickness.  Figure  11  shows  the  MTF  data  in  Fig.  10(c)  plot¬ 
ted  as  a  function  of  the  dimensionless  quantity  fL,  where 
/  is  the  spatial  frequency  and  L  is  the  physical  thickness  of 
the  Csl  layers.  In  order  to  reflect  the  actual  mass  loading 
of  different  Csl  layers,  the  estimated  thickness  shown  in 
Table  1  was  used  instead  of  the  nominal  thickness.  The 
universal  MTF  plot  was  first  introduced  by  Swank  in  order 
to  show  that  when  plotted  this  way,  the  MTF  is  the  same  for 
all  powder  phosphor  screens  of  the  same  design.10  Also  plot¬ 
ted  in  Fig.  1 1  is  the  universal  MTF  of  a  powder  phosphor 
screen  (Kodak  Lanex  Regular)  which  has  a  thickness  of  1 67 
fi m.1  It  shows  that  the  MTF  for  Csl  of  the  same  thickness  is 
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Fig.  10.  (a)  Presampling  MTF  of  the  Csl  layers  for  all  the  Csl  samples;  (b) 
Theoretical  calculation  of  MTFK,  i.e.,  MTF  due  to  K-fluorescence  reabsorp¬ 
tion  for  each  Csl  layer;  and  (c)  estimated  MTF0  of  the  Csl  layers,  i.e.,  MTF 
due  to  optical  factors,  by  dividing  the  MTF  data  in  Fig.  10(a)  by  the  calcu¬ 
lated  MTFk  shown  in  Fig.  10(b). 


more  than  twice  as  high  as  the  powder  phosphor,  demonstrat¬ 
ing  the  advantage  of  the  columnar  structure.  Figure  1 1  also 
shows  that  the  universal  MTF  for  Csl  increases  as  a  function 
of  L,  indicating  that  the  benefit  of  the  columnar  structure  is 


fL 


Fig.  11.  Universal  MTF  plots  for  the  HL  Csl  samples,  i.e.,  MTF  plotted  as 
a  function  of  the  product  of  the  spatial  frequency  / and  the  thickness  of  Csl 
L.  Since  the  optical  characteristics  of  the  HR  type  is  different  from  that  of 
the  HL,  only  the  MTF  for  HL  Csl  layers  are  plotted  to  emphasize  the 
thickness  dependence.  Plotted  in  the  same  graph  for  comparison  is  the  uni¬ 
versal  MTF  plot  for  a  Kodak  Lanex  Regular  powder  phosphor  screen  (re¬ 
plotted  from  Ref.  1). 

relatively  more  pronounced  in  thicker  Csl  layers.  This  is  con¬ 
sistent  with  our  hypothesis  that  each  Csl  sample  consists  of  a 
structureless  layer  with  a  fixed  thickness,  which  makes  the 
fraction  of  columnar  structure  increases  as  a  function  of 
thickness. 

C.  DQE 

In  order  to  determine  the  DQE  experimentally,  the  x-ray 
sensitivity  and  NPS  measurements  at  different  exposures  are 
needed  in  addition  to  the  MTF.  Shown  in  Fig.  12  is  the 
measured  x-ray  sensitivity  of  the  150  /xm  HL  Csl  layers  with 
x-ray  exposures  from  0.14  to  20.9  mR.  The  same  graph 
shows  the  best  linear  fit  of  the  experimental  data.  Nonlinear¬ 
ity  is  observed  in  the  detector  response.  Similar  results  were 
obtained  from  other  Csl  samples.  We  believe  the  nonlinearity 
is  due  to  the  response  of  the  CMOS  optical  sensor,  which  has 
been  reported  previously,17  although  not  as  severe  as  ob¬ 
served  in  our  experiments.  The  following  approach  was 
taken  to  correct  for  the  effect  of  nonlinearity  in  the  measure¬ 
ment  of  NPS:  (1)  performed  flat-fielding  of  noise  images  as 
described  in  Sec.  II C;  (2)  at  each  x-ray  exposure,  performed 
best  linear  fitting  of  the  x-ray  sensitivity  with  3  data  points 
(including  the  exposures  just  below  and  just  above)  using 
function  y^a  +  bx;  (3)  corrected  the  flat-fielded  images  us¬ 
ing: 

FFj  j— a 

CIiJ=—Ji - ko,  (5) 

where  C/;J  and  FFtj  were  the  nonlinearity  corrected  noise 
image  and  flat-fielded  image  with  pixel  coordinate  of  i  and  j, 
respectively,  and  k0  was  the  overall  x-ray  sensitivity  for  the 
entire  exposure  range  obtained  from  the  best  linear  fitting  as 
shown  in  Fig.  12;  (4)  used  Fourier  transform  to  calculate  2D 
NPS  as  described  in  Sec.  II C.  Figure  13  shows  the  calcu- 
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Fig.  12.  An  example  of  x-ray  linearity  measurements:  Pixel  response  plotted 
as  a  function  of  x-ray  exposure  for  the  150  /xm  HL  Csl  layer. 


lated  NPS  at  four  exposure  levels  for  the  150  /xm  HR  Csl 
layer,  with  0.26  mR  equivalent  to  the  mean  exposure  to  the 
detector  in  general  radiography.  Also  shown  is  the  measure¬ 
ment  of  the  electronic  noise  power  of  the  system  (due  to  the 
electronic  noise  of  the  CMOS  optical  sensor).  It  shows  that 
the  electronic  noise  contributes  significantly  to  the  NPS  at 
high  spatial  frequencies.  Shown  in  Fig.  14(a)  are  the  calcu¬ 
lated  DQE(/)  of  the  same  Csl  layer.  The  results  showed  that 
DQE(0)  measured  at  all  exposures  are  essentially  identical. 
However  at  high  spatial  frequencies,  DQE  is  compromised 
by  the  electronic  noise,  especially  at  low  exposures,  which 
indicates  that  the  low  conversion  of  HR  type  Csl  makes  it 
more  susceptible  to  electronic  noise.  For  comparison,  in  Fig. 
14(b)  is  the  DQE  measurement  for  the  150  /xm  HL  Csl  layer, 
which  has  higher  conversion  gain  than  the  HR  type.  It  shows 
that  the  DQE  at  high  spatial  frequencies  reaches  a  limit  for 
exposures  of  1 .08  mR  and  larger.  The  DQE(0)  of  both  the 
HR  and  HL  150  /xm  layers  are  —0.37,  which  confirms  that 
the  150  /xm  HR  and  HL  layers  have  similar  x-ray  absorption 
(mass  loading)  and  Swank  factor.  These  extrapolated 
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Fig.  13.  NPS  measurements  in  the  dark  and  at  four  different  x-ray  exposure 
levels  for  the  150  /xm  HR  Csl  layer. 
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Fig.  14.  DQE  measurements  at  four  different  x-ray  exposures  for  the  150 
/xm  Csl  layers:  (a)  HR;  and  (b)  HL  types. 


DQE(0)  values  from  the  DQE(/)  measurements  also  agree 
well  with  the  independent  PHS  measurement  of  DQE(0) 
shown  in  Fig.  9  for  the  mean  energy  of  54  keV  of  the  spec¬ 
trum  used  in  our  imaging  experiments. 

In  order  to  determine  the  theoretical  limit  of  the  DQE  of 
different  Csl  layers  used  in  our  experiments,  and  facilitate 
the  comparison  with  DQE  models,  we  subtracted  the  elec¬ 
tronic  noise  from  the  NPS  measurements  at  different  expo¬ 
sures.  This  makes  the  DQE  for  each  Csl  layer  exposure  in¬ 
dependent  (within  experimental  error).  For  clarity  only  the 
subtracted  results  for  the  highest  exposure  are  shown  in 
Fig.  15. 

D.  DQE  model 

The  results  of  DQE  calculation  using  the  linear  system 
model  are  shown  in  Fig.  15  along  with  experimental  mea¬ 
surements  for  comparison.  The  measured  MTF  and  mass 
loading  of  each  Csl  layer  were  used  for  the  calculation  of 
DQE,  and  the  model  did  not  include  the  Lubberts  effect.  As 
shown  in  Fig.  15  for  the  150  /xm  HR  Csl,  the  modeled  and 
measured  DQE  agree  very  well,  which  means  that  there  is 
very  little  Lubberts  effect  in  HR  type  Csl  layer.  For  the  HL 
Csl  layers,  the  measured  and  modeled  DQE  agree  well  at 
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Spatial  frequency  (cycles/mm)  Spatial  frequency  (cycles/mm) 

Fig.  15.  Comparison  between  the  measured  and  modeled  DQE  for  four  different  Csl  layers:  The  150  fim  HR  and  HL  layers,  the  300  fim  and  the  600  fim  HL 
layers.  (Since  the  mass  loading  and  MTF  for  the  450  and  600  fim  HL  Csl  layers  are  very  similar,  only  the  result  for  the  600  ^m  layer  is  shown.)  After  the 
electronic  noise  is  subtracted  from  the  measured  DQE  at  different  exposures,  the  DQE  is  no  longer  exposure  dependent  (within  experimental  errors).  Hence 
for  clarity  only  the  subtracted  results  for  the  highest  exposure  level  are  shown. 


low  spatial  frequencies.  However  at  high  spatial  frequencies, 
the  modeled  DQE  is  considerably  higher  than  the  measured 
results.  This  indicates  Lubberts  effect  in  the  HL  type  Csl. 
The  difference  between  measured  and  modeled  results  in¬ 
creases  as  a  function  of  the  thickness  of  Csl,  suggesting  Lub¬ 
berts  effect  is  more  severe  in  thicker  HL  Csl  layers.  Since  the 
Lubberts  effect  arises  as  a  consequence  of  different  MTF  for 
x-rays  absorbed  at  different  depths  within  the  Csl,  it  can  be 
concluded  that  these  MTF  variations  are  lowest  for  150  yum 
HR  layer  and  then  increase  both  as  the  HL  structure  is 
adopted  and  then  even  further  as  layer  thickness  is  increased. 

Our  results  of  DQE  measurement  and  modeling  indicate 
that  compared  to  the  widely  used  HL  Csl,  HR  Csl  has  less 
Lubberts  effect,  which  can  lead  to  improved  DQE  at  high 
spatial  frequencies.  However  this  is  at  the  cost  of  greater 
susceptibility  to  electronic  noise  and  compromised  imaging 
performance  at  low  dose.  If  additional  gain  can  be  obtained 
in  indirect  flat-panel  detectors,  e.g.,  through  avalanche  gain, 
HR  Csl  layers  could  be  used  to  boost  the  DQE  at  high  spatial 
frequencies.28  This  is  particularly  important  for  high  resolu¬ 
tion  x-ray  imaging  applications,  e.g.,  mammography. 

Since  our  experiments  are  based  on  Csl  samples  from  one 


manufacturer,  who  keeps  many  details  of  how  the  samples 
are  prepared  proprietary,  it  may  be  difficult  to  apply  our  find¬ 
ings  directly  to  Csl  made  by  other  manufactures.  Our  hope  is 
that  through  the  establishment  of  a  standard  methodology  for 
investigating  Csl  imaging  properties,  similar  measurements 
can  be  performed  by  others  for  comparison.  We  also  hope 
that  our  understanding  of  the  physics  behind  the  dependence 
of  imaging  properties  on  thickness  and  design  types  can  lead 
to  better  modeling  of  the  light  channeling  process  in  Csl,  and 
provide  a  source  of  experimental  data  for  model  validation. 

IV.  CONCLUSIONS 

In  this  paper,  we  presented  our  experimental  study  of  the 
imaging  performance  of  structured  Csl.  Both  HR  and  HL 
type  Csl,  as  well  as  HL  Csl  of  different  thicknesses  Were 
investigated.  Our  results  showed  that  HR  and  HL  Csl  layers 
of  the  same  thickness  have  essentially  the  same  Swank  fac¬ 
tor.  The  HR  Csl  has  less  Lubberts  effect  compared  to  the  HL, 
but  at  the  cost  of  lower  conversion  gain,  which  leads  to  a 
greater  susceptibility  to  electronic  noise.  The  MTF  of  all  the 
samples  dropped  rapidly  at  low  spatial  frequencies,  and  with 
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a  magnitude  which  was  identical  for  all  thicknesses.  This 
suggests  a  similar  components  in  all  the  samples  which  was 
hypothesized  to  be  a  transparent  (unstructured)  layer  close  to 
the  FOP  substrate.  Although  MTF  of  the  HL  Csl  decreases  as 
a  function  of  thickness,  on  the  universal  MTF  plot,  i.e.,  MTF 
plotted  as  a  function  of  the  product  of  spatial  frequency  / and 
thickness  L,  MTF  increases  with  thickness.  This  demon¬ 
strates  that  the  advantage  of  columnar  structure  in  Csl  is 
more  pronounced  for  thicker  layers. 
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ABSTRACT 

An  indirect  flat-panel  imager  (FPI)  with  avalanche  gain  is  being  investigated  for  low-dose  x-ray  imaging.  It  is  made  by 
optically  coupling  a  structured  x-ray  scintillator  CsI(Tl)  to  an  amorphous  selenium  (a-Se)  avalanche  photoconductor 
called  HARP.  The  final  electronic  image  can  be  read  out  using  either  an  array  of  thin  film  transistors  (TFT)  or  field 
emitters  (FE).  The  advantage  of  the  proposed  detector  is  its  programmable  gain,  which  can  be  turned  on  during  low  dose 
fluoroscopy  to  overcome  electronic  noise,  and  turned  off  during  high  dose  radiography  to  avoid  pixel  saturation.  This 
paper  investigates  the  important  design  considerations  for  HARP  such  as  avalanche  gain,  which  depends  on  both  the 
thickness  dsc  and  the  applied  electric  field  E$e.  To  determine  the  optimal  design  parameter  and  operational  conditions  for 
HARP,  we  measured  the  ESe  dependence  of  both  avalanche  gain  and  optical  quantum  efficiency  of  an  8  pm  HARP  layer. 
The  results  were  applied  to  a  physical  model  of  HARP  as  well  as  a  linear  cascaded  model  of  the  FPI  to  determine  the 
following  x-ray  imaging  properties  in  both  the  avalanche  and  non-avalanche  modes  as  a  function  oiESe'  (1)  total  gain 
(which  is  the  product  of  avalanche  gain  and  optical  quantum  efficiency);  (2)  linearity;  (3)  dynamic  range;  and  (4)  gain 
non-uniformity  resulting  from  thickness  non-uniformity.  Our  results  showed  that  a  HARP  layer  thickness  of  8  pm  can 
provide  adequate  avalanche  gain  and  sufficient  dynamic  range  for  x-ray  imaging  applications  to  permit  quantum  limited 
operation  over  the  range  of  exposures  needed  for  radiography  and  fluoroscopy. 

1.  INTRODUCTION 

We  previously  proposed  the  concept  of  indirect  flat-panel  imagers  (FPI)  with  avalanche  gain  by  optically  coupling  a 
structured  x-ray  scintillator  CsI(Tl)  to  an  amorphous  selenium  (a-Se)  avalanche  photoconductor  called  HARP  (High-gain 
Avalanche  Rushing  Photoconductor). 1  Shown  in  Figure  1  is  a  side  view  of  the  proposed  detector.  The  photons  emitted 
from  the  Csl  are  absorbed  by  the  HARP  layer  and  generate  electron-hole  pairs  near  the  top  interface  (light  entrance  side). 
Under  a  sufficiently  high  electric  field  ESe  (i.e.  applying  a  positive  voltage  on  the  top  bias  electrode  of  a-Se),  holes  move 
towards  the  bottom  surface  and  undergo  avalanche  multiplication.  The  avalanche  gain  gav  can  be  as  high  as  1500 
depending  on  Ese  and  the  thickness  of  the  a-Se  layer,  dSe ■  The  holes  are  collected  at  the  bottom  surface  of  a-Se  to  form  a 
charge  image,  which  is  read  out  electronically  using  either  an  active  matrix  thin  film  transistor  (TFT)  array  or  a  field 
emitters  array  (FEA).  We  called  the  proposed  detector  with  TFT  readout  SHARP- AMFPI  (Scintillator  HARP  Active 
Matrix  Flat  Panel  Imager)  and  with  FEA  readout  S  APHIRE  (Scintillator  Avalanche  Photoconductor  with  High 
resolution  Emitter  Readout).  The  advantage  of  the  proposed  detectors  is  its  programmable  gain,  which  can  be  turned  on 
during  low  dose  imaging  procedures  (e.g.  fluoroscopy)  to  overcome  electronic  noise,  and  turned  off  during  high  dose 
imaging  (e.g.  radiography)  to  avoid  pixel  saturation.  This  paper  investigates  the  important  design  parameters  and 
operational  conditions  for  HARP  in  order  to  ensure  adequate  dynamic  range  for  both  fluoroscopy  and  radiography. 
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Figure  1 :  Side  view  showing  the  concept  of  indirect  flat-panel  detector  with  avalanche  gain.  (Note:  thickness  not  to  scale.) 
X-rays  enter  from  the  top  with  the  ITO  (top)  electrode  of  HARP  biased  positively.  The  charge  (hole)  image  collected  on  the 
bottom  surface  of  the  HARP  is  read  out  with  either  a  TFT  array  or  a  FEA. 

Since  avalanche  multiplication  in  a-Se  was  discovered  over  a  decade  ago,  HARP  tubes,  which  consist  of  a  HARP  layer 
read  out  with  a  scanning  electron  beam,  have  been  commercialized  for  high  sensitivity  and  high  definition  television 
applications.  Stable  and  uniform  avalanche  multiplication  has  been  observed  with  very  little  added  noise.2  The  electric 
field  ESe  under  which  avalanche  multiplication  occurs  is  >  70  V/pm,  which  is  an  order  of  magnitude  higher  than  that 
typically  used  in  direct  FPI  incorporating  thick  a-Se  layers.  To  sustain  an  ESe  of  this  magnitude  without  a  significant 
increase  in  dark  current  due  to  charge  injection  from  the  bias  electrodes,  blocking  layers  need  to  be  developed  for  both 
electrode  interfaces.  As  shown  in  Figure  1,  the  HARP  sensor  consists  of  a  multi-layer  structure.3  It  contains  a  transparent 
bias  electrode  indium  tin  oxide  (ITO)  that  is  biased  positively  during  imaging,  followed  by  a  thin  layer  (~20  nm)  of 
Cerium  Oxide  (Ce02)  and  a  thin  layer  of  lithium  fluoride  (LiF)  doped  a-Se  to  block  the  injection  of  holes  from  the  ITO 
bias  electrode.  The  intrinsic  a-Se  photoconductor  is  usually  0.5  -  35  pm  thick  depending  on  the  desired  avalanche  gain 
and  doped  with  arsenic  (As)  to  prevent  the  a-Se  from  recrystallizing.  The  bottom  interface  of  the  HARP  structure  is  a 
sub-micron  thin  layer  of  antimony  tri-sulfide  (Sb2S3),  which  is  used  to  block  electron  injection  from  the  negative  bias 
electrode  (in  a  solid  state  detector  with  TFT  readout)  or  the  scanning  electron  beam  (in  FIARP  tubes  or  with  FEA 
readout).  Since  the  mobility  of  holes  in  a-Se  is  more  than  an  order  of  magnitude  higher  than  that  of  the  electrons,  only 
holes  create  a  significant  number  of  new  carriers  by  impact  ionization.  Hence  it  is  essential  for  light  photons  to  enter 
from  the  positive  bias  side. 

Our  previous  calculation  using  a  preliminary  linear  system  model  of  indirect  FPI  with  avalanche  gain  predicted  that  a 
gain  of  20  would  be  needed  to  produce  x-ray  quantum  noise  limited  performance  for  fluoroscopy  and  low  dose 
mammography  (tomosynthesis).  This  requires  a  HARP  layer  thickness  of  8  pm  or  more.  To  determine  proper  design 
parameters  and  operational  conditions  for  HARP,  we  measured  the  ESe  dependence  of  both  avalanche  gain  and  optical 
quantum  efficiency  of  an  8  pm  HARP  layer.  The  results  were  used  in  a  physical  model  of  HARP  as  well  as  a  linear 
cascaded  model  of  the  FPI  to  determine  the  following  x-ray  imaging  properties  in  both  the  avalanche  and  non-avalanche 
modes  as  a  function  of  dSe  and  ESe:  (1)  total  gain  (which  is  the  product  of  avalanche  gain  and  optical  quantum  efficiency); 
(2)  linearity;  (3)  dynamic  range;  and  (4)  gain  non-uniformity  resulting  from  thickness  non-uniformity. 

2.  PHOTOSENSITIVTY  OF  HARP 

The  x-ray  imaging  performance  of  indirect  FPI  with  avalanche  gain  relies  on  the  photosensitivity  of  HARP. 
Understanding  its  dependence  on  ESe  and  layer  thickness  dSe  is  crucial  for  optimization  of  HARP  design  parameters  and 
operational  conditions  for  both  the  avalanche  and  non-avalanche  modes.  In  this  section,  we  will  first  present  the 
experimental  method  and  results  of  the  photosensitivity  measurements  of  HARP,  and  then  describe  the  theory  of  ESe 
dependence  of  photosensitivity,  which  is  proportional  to  both  the  optical  quantum  efficiency  (QE)  q  and  the  avalanche 
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gain  gm.  The  theory  and  experimental  data  will  then  be  used  to  determine  r|  and  gav  as  a  function  of  ESe  in  the  context  of 
x-ray  imaging. 

2.1.  Experimental  measurement  of  photosensitivity 

2.1.1.  Method: 


The  experimental  setup  for  the  photosensitivity  measurement  is  shown  in  Figure  2.  A  grating  monochromator  was  used 
to  generate  optical  photons  of  a  single  wavelength  X.  The  range  of  X  values  used  in  our  experiments  was  between  300- 
700  nm.  The  intensity  of  the  output  beam  was  attenuated  by  neutral-density  (ND)  filters.  A  beamsplitter  was  used  to 
direct  a  fraction  (90  %)  of  the  attenuated  beam  to  a  silicon  (Si)  photodiode  (PD),  which  monitored  the  beam  intensity 
during  experiment.  The  optical  photons  passing  directly  through  the  beam  splitter  (10  %)  were  detected  by  the  HARP 
tube.  Before  sensitivity  measurements  the  HARP  input  light  intensity  was  calibrated  for  each  X  in  order  to  keep  a 
constant  input  power.  This  was  achieved  by  adjusting  the  ND  filters  and  monitoring  the  light  intensity  using  the  Si  PD. 


For  each  X  and  ESe,  the  signal  current  S  obtained  by  the  HARP  tube  was  recorded,  and  the  photosensitivity  was 
quantified  as  the  number  of  signal  charge  carriers  generated  by  one  absorbed  light  photon.  This  quantity  is  referred  to  as 
the  effective  quantum  efficiency  r)  ,  and  it  is  the  product  of  q  and  gav.  The  values  of  T|*  was  experimentally  determined 
from  the  measurement  of  the  HARP  tube  current  S  using: 


S/e 

IT/hv 


(2) 


where  e  was  the  electron  charge,  /  was  the  HARP  input  light  intensity  power  given  in  Watts  (W),  hv was  the  energy  of 
each  incident  photon,  and  T  was  the  fraction  of  the  incident  photons  (to  HARP)  reaching  the  a-Se  layer.  The  correction 
factor  T  was  due  to  light  attenuation  by  both  the  ITO  bias  electrode  and  the  Ce02  blocking  layer  for  HARP,  and  is  a 
function  of  the  wavelength  X.  The  value  of  T  was  experimentally  measured  for  all  X  used. 


2.1.2.  Measurement  results: 

Shown  in  Figure  3  is  the  measured  effective  quantum  efficiency  q*  as  a  function  of  ESe  for  an  8  pm  thick  HARP  layer.  It 
shows  that  for  blue  light  (X=400  nm),  q  reaches  a  plateau  at  Ese  >  20  V/pm.  This  corresponds  to  a  saturation  of  optical 
quantum  efficiency  q  at  unity.  With  further  increase  in  ESe  to  >  80  V/pm,  q  starts  to  increase  again,  demonstrating  the 
onset  of  avalanche.  However  for  green  light,  the  ESe  dependence  of  q*  is  quite  different.  As  shown  in  Figure  3  for  X=540 
nm,  which  is  the  peak  of  the  emission  spectrum  of  Csl  (Tl),  q  begins  at  a  much  smaller  value  compared  to  X=400  nm.  It 
then  increases  continuously  with  ESe  without  saturation  before  avalanche  starts  at  ~80  V/pm.  This  is  because  of  the  lower 
optical  q  for  X=540  nm  compared  to  X=400  nm.  It  does  not  saturate  as  a  function  of  ESe  or  reach  unity  until  ESe  is  well 
into  the  range  of  avalanche  multiplication.  To  choose  the  proper  operational  condition  for  HARP  in  different  x-ray 
imaging  applications,  we  need  to  understand  separately  the  ESe  dependence  of  q  and  gm,  This  will  be  discussed  in  section 

2.2. 
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Figure  3:  Measured  effective  quantum  efficiency  t|*  as  a  function  of  ESc  for  an  8  pm  thick  HARP  layer. 

2.2.  Theory  of  photosensitivity  of  HARP:  Quantum  efficiency  and  avalanche  gain 


Figure  4:  The  optical  quantum  efficiency  of  a-Se  calculated  using  Onsager  theory  and  r0-  1.7  nm,  compared  with  the 
measured  effective  quantum  efficiency  ri*  at  X  =  540  nm. 


The  effective  quantum  efficiency  13*  of  HARP  is  proportional  to  both  13  and  gav.  The  value  of  13  in  a-Se  depends  on  both 
ESe  and  the  wavelength  X  of  the  incident  photon.  It  was  found  by  Pai  and  Enck  that  the  ESe  and  X  dependence  of 
photogeneration  in  a-Se  follows  the  Onsager  mechanism.4  Onsager  theory  states  that  every  absorbed  photon  creates  a 
pair  of  thermalized  carriers  bound  by  their  Coulombic  attraction.  The  initial  separation  r0  of  the  charge  pair  increases 
with  incident  photon  energy.  The  probability  of  dissociation  for  the  bound  pairs  into  free  carriers  increases  with  r0  and 
the  applied  external  field  ESe  with  the  relationship: 


,o  kT  c-*c-'E*n>i*r  y 

0  eESer0  ml 


m= 0  ' 


y  y  I,  with 

tSz-Sif.  kT  ’  U  4 nee.kT 


(1) 
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where  <j>0  is  defined  as  the  efficiency  of  production  of  thermalized  charge  pairs  per  absorbed  photon,  and  it  is  independent 
of  Ese-  Pai  and  Enck  showed  that  the  room  temperature  measurement  of  r)  can  be  best  fitted  with  the  expression  in  Eq.  1 
using  <(>o=l  and  r0  in  the  range  between  0.8  -  7.0  nm  for  A.  values  between  620  and  400  nm.  From  now  on  we  will  focus 
on  A=540  nm,  which  is  the  peak  of  the  emission  spectrum  for  Csl  (Tl).  Shown  in  Figure  4  is  the  Tj  calculated  using  Eq.  1 
and  r0  =  1 .7  nm.  Also  plotted  in  Figure  4  is  the  measured  r\  at  the  same  wavelength  for  dSe  —  8  jam.  Figure  4  shows  that 
the  Onsager  theory  provides  a  reasonable  prediction  for  the  ESe  dependence  of  T|  in  HARP  for  ESe  <  70  where  gav  =  1 . 
The  calculated  t|  using  Onsager  theory  was  then  divided  into  the  measured  Tj*  so  that  the  £&  dependence  of  avalanche 
gain  gav  could  be  obtained.  The  results  are  shown  in  Figure  5  (a),  where  gav  =  46  at  ESe  =110  V/pm.  Figure  5  (a)  shows 
that  gav  increases  very  rapidly  with  ESe.  Hence  choosing  the  right  operating  condition  for  HARP  is  crucial.  To  determine 
the  x-ray  response  of  the  indirect  FPI  with  avalanche  gain,  it  is  necessary  to  understand  the  continuous  dependence  of  gav 
on  ESe.  This  is  because  as  image  charge  accumulates  on  the  pixels,  the  voltage  drop  across  the  HARP  layer  will  decrease. 
With  the  rapid  change  in  gav  as  a  function  of  ESe  shown  in  Figure  5  (a),  the  value  of  gaY  will  decrease  as  a  result  of  charge 
accumulation  on  the  pixels.  This  will  affect  the  x-ray  sensitivity,  linearity  and  dynamic  range  of  the  FPI. 


e^fV/nm)  l/E^m/V) 

(a)  (b) 


Figure  5:  (a)  Solid  circles  are  the  gav  for  the  8  pm  HARP  layer  calculated  by  dividing  the  measured  Tj*  at  540  nm  by  the 
optical  Tj  predicted  by  the  Onsager  theory  using  r0  =  1.7  nm.  The  solid  line  shows  the  gav  calculated  using  the  fitted  pi  and  p2 
values,  (b)  The  plot  of  ln(ln(g„v))  as  a  function  of  1  /ESe,  and  the  best  linear  fit  for  the  data. 


The  avalanche  gain  characteristics  in  HARP  have  been  investigated  as  a  function  of  dSe  and  ESe.5  According  to  the  theory 
of  avalanche  multiplication,6  when  ESe  across  the  thickness  dSe  is  uniform,  and  the  distance  through  which  avalanche 
occurs  is  equal  to  dSe  (which  is  true  for  HARP  because  photons  are  absorbed  within  0.1  pm  at  the  light  entrance  side  of 
the  a-Se  layer),  gav  is  given  by: 

=  (/?-g0^ 

8m  p-cced^-a) 

where  a  and  p  are  the  impact  ionization  rates  of  electrons  and  holes,  respectively.  Since  in  a-Se  the  drift  mobility  of 
holes  (0.18  cm2/ Vs)  is  much  higher  than  that  of  electrons  (0.003  cm2/Vs),  p  is  much  greater  than  a.  Thus  in  HARP  light 
always  enters  from  the  positive  bias  electrode  to  maximize  the  avalanche  gain.  Under  the  assumption  that  the  impact 
ionization  of  electrons  is  negligible,  Eq.  2  can  be  simplified  to: 

gm=e“^  (3) 


with  p  given  by: 


dh. 


(4) 


The  values  for  P]  and  p2  can  be  determined  by  curve  fitting  of  the  experimental  measurement  of  gav.  By  taking  the 
logarithm  of  Eq.  3  twice,  we  obtain: 


356  Proc.  ofSPIEVol.  5745 


(5) 


ln(ln  gav)  =  \n(PxdSe)  , 

^Se 

which  indicates  that  the  quantity  ln(lngav)  is  linearly  related  to  1  /ESe.  Shown  in  Figure  5  (b)  is  the  experimentally 
determined  gav  data  [shown  in  Figure  5(a)]  replotted  in  the  form  of  ln(lngm,)  as  a  function  of  l/ESe.  The  best  linear  fit  for 
the  experimental  data  is  also  shown  in  Figure  5  (b).  From  the  intercept  and  the  slope  of  the  fitted  curve  we  obtained 
pi=5.5  x  103  /pm  and  p2  =  1.029  x  103  V/pm.  Using  these  pt  and  p2  values  and  the  analytical  expression  for  gav  shown  in 
Eqs.  3  and  4,  the  theoretical  gov  values  are  calculated  and  shown  in  Figure  5  (a).  The  calculated  gav  using  the  fitted 
parameters  agrees  well  with  experimentally  determined  gav  except  for  the  highest  ESe  of  1 12.5  V/pm.  This  could  be  seen 
from  the  relative  large  error  in  fitting  for  this  data  point  in  Figure  5  (b).  Since  gm,  is  related  to  p2  through  two  exponential 
functions,  a  small  error  in  the  estimate  of  p2  can  result  in  a  large  error  in  the  estimate  of  gav.  Nevertheless  the  quality  of 
fitting  for  ESe  2  110  V/pm  is  excellent,  and  the  fitted  parameters  pi  and  p2  will  be  used  in  the  following  calculation  of  x- 
ray  response. 


3.  X-RAY  IMAGING  PERFORMANCE 


We  have  previously  developed  a  cascaded  linear  system  model  for  the  indirect  FPI  with  avalanche  gain.1  In  this  section, 
we  will  apply  the  ESe  dependence  of  gov  and  q  derived  in  the  previous  section  to  the  cascaded  linear  system  model  to 
determine  the  x-ray  response  and  the  detective  quantum  efficiency,  DQE(/),  of  the  indirect  FPI  with  avalanche  gain.  We 
will  use  the  resolution  and  avalanche  gain  requirements  for  a  radiography/fluoroscopy  (R/F)  detector  as  an  example.  It 
can  be  shown  from  the  cascaded  linear  system  model  that  the  mean  signal  from  each  pixel  of  the  detector  is  given  by:1 

<Ps=ap2a7igavYJTlxg(E)qQ(E)  (6) 

E 

where  ap  is  the  pixel  pitch  of  the  detector,  qr  is  the  x-ray  quantum  absorption  efficiency  of  the  Csl  layer,  g(E)  is  the  x-ray 
to  light  conversion  gain  of  Csl  that  depends  on  the  x-ray  energy  E  and  the  energy  required  to  generate  an  optical  photon 
W,  qo(E)  is  the  incident  x-ray  spectrum  (which  is  given  in  number  of  photons  per  unit  area  for  a  given  exposure),  and  a  is 
the  optical  coupling  efficiency  between  the  Csl  and  the  a-Se  layer.  The  DQE(/)  of  the  detector  is  given  by:1 

T'a  (f)Tb  (f)Tc2  if) 


DQE  (/)  =  - 


[~~T2(f)T2(f)  + 
?hAs 


^gav  ^)  -1'T’  2 


(7) 


gmgVx<xn 


F/CO+ 


(ap2gmgnxO(T])2q(s 


where  As  is  the  Swank  factor  of  Csl,  Sa  is  the  added  electronic  noise,  TJf),  Tb(f)  and  Tc(f)  are  the  modulation  transfer 
function  (MTF)  due  to  the  Csl,  the  optical  coupling  between  Csl  and  HARP,  and  the  pixel  aperture  function,  respectively. 


Our  previous  investigation  of  the  inherent  x-ray  imaging  performance  of  Csl  (Tl)  has  shown  that  a  600  pm  Csl  layer 
optimized  for  high  light  output  (HL)  could  provide  adequate  performance  for  R/F  applications.7  Hence  the  measured 
imaging  properties  of  this  Csl  layer  were  used  in  our  cascaded  linear  system  model.  Table  I  summarizes  the  detector 
parameters  and  HARP  operating  conditions  chosen  for  fluoroscopy  and  radiography. 


Table  I:  Detector  parameter  and  operating  conditions  chosen  for  different  x-ray  imaging  applications. 


R/F  detector 

Fluoroscopy 

Radiography  | 

Pixel  pitch  aD  (pm) 

200 

X-ray  spectrum 

RQA5 

Electronic  noise 

1500  electrons 

Csl  (75%  packing  density) 

600  pm  HL 

If  for  Csl  (eV) 

18 

a 

0.8 

Minimum  exposure 

0.1  pR 

30  pR 

ESe  (V/pm) 

110 

95 

A 

0.36 

0.31 

gav 

46 

2.3 
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3.1.  DQE (f) 


Shown  in  Figure  6  (a)  is  the  DQE(/)  calculated  using  the  cascaded  linear  system  model  for  the  detector  with  different  gm, 
settings.  The  calculation  was  performed  for  the  lowest  exposure  in  fluoroscopy  (0. 1  pR),  where  the  detector  performance 
is  most  susceptible  to  electronic  noise.  It  shows  that  gm,  =  46  is  sufficient  to  bring  the  DQE(/)  to  its  theoretical  limit  (i.e.  a 
detector  without  electronic  noise).  The  DQE(/)  for  the  next  highest  gav  =  1 2  measured  at  ESe  =  1 05  V/gm  also  provides 
reasonable  DQE.  Thus  an  ESe  between  105  and  1 10  V/|im  will  be  adequate  for  fluoroscopy.  Figure  6  (b)  shows  the 
DQE(/)  calculated  for  the  lowest  radiographic  exposure  of  30  pR.  It  shows  that  ESe  =  95  V/pm  can  provide  sufficient  p 
and  gav.  Hence  the  ESe  chosen  in  Table  I  is  adequate  for  the  dual-mode  operation  of  an  RTF  detector. 


Spatial  frequency  (cycles/mm) 
(a) 


Spatial  frequency  (cycles/mm) 
(b) 


Figure  6:  Calculated  DQE(/)  using  the  detector  parameters  and  operating  conditions  shown  in  Table  I  for  a  fluoroscopy 
detector  at  the  x-ray  exposure  level  of:  (a)  0.1  pR;  and  (b)  30  pR. 

3.2.  Pixel  x-ray  response 

Shown  in  Figure  7  is  the  calculated  pixel  x-ray  response  for  the  R/F  detector  under  two  different  operational  conditions. 
It  was  assumed  that  the  dominant  pixel  capacitance  was  from  the  capacitance  of  the  HARP  layer.  This  condition  results 
in  the  greatest  drop  in  the  potential  across  the  HARP  layer  for  a  given  x-ray  exposure,  and  creates  the  largest  effect  on 
linearity.  Figure  7  (a)  shows  that  for  fluoroscopy,  the  initial  ESe=  110  V/pm  provides  a  gm,  of  46.  Within  the  regular 
exposure  range  of  0. 1  -  10  pR/frame  in  fluoroscopy,  the  detector  has  very  good  linearity,  i.e.  an  essentially  constant  gav. 
With  further  increase  in  exposure,  which  corresponds  to  detector  exposure  around  the  edge  of  the  human  body  or  raw 
exposure  to  detector,  the  x-ray  response  becomes  sublinear  with  a  decrease  in  gav.  This  decrease  in  gav  ensures  a  wide 
dynamic  range  for  the  detector  because  there  is  no  “hard”  saturation.  The  pixel  potential  will  continue  to  increase  as  a 
function  of  x-ray  exposure  until  it  reaches  the  predetermined  Vmax  for  high  voltage  protection  (for  the  TFT  readout 
method).  The  method  for  high  voltage  protection  is  described  in  detail  elsewhere8  and  will  not  be  discussed  here.  The 
pixel  potential  for  the  exposure  of  1  mR  in  the  fluoroscopy  mode  is  VP  =  56  V,  which  is  safe  for  a-Si  TFTs.  Figure  7  (b) 
shows  that  with  a  much  lower  gav  =  2.4  programmed  for  radiography,  the  detector  response  for  the  regular  detector 
exposure  range  of  30  -  3  mR  is  essentially  linear  with  very  little  change  in  gav.  Comparing  Figures  7  (a)  and  (b),  the 
image  charge  on  the  pixel  is  in  the  same  range  despite  their  large  difference  in  exposure  range.  This  demonstrates  the 
wide  dynamic  range  (over  5  orders  of  magnitude)  of  the  detector  by  varying  gav  for  different  x-ray  imaging  applications. 
This  is  different  from  a  detector  with  constant  high  gain,  which  has  the  problem  with  pixel  saturation  for  radiographic 
applications. 
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Figure  7:  The  calculated  image  charge  on  each  pixel  of  the  detector  (left  axis)  and  the  corresponding  avalanche  gain  (right 
axis)  as  a  function  of  x-ray  exposure  for  the  detector  operational  setting  of:  (a)  fluoroscopy;  and  (b)  radiography 

4.  EFFECT  OF  HARP  THICKNESS  UNIFORMITY 


An  important  issue  for  making  a  large  area  HARP  is  its  thickness  uniformity.  Currently  HARP  films  can  be  made  very 
uniform  for  an  area  of  1-2”  in  diameter  with  essentially  no  visible  gain  variation  in  the  images.  The  images  currently 
produced  by  HARP  for  broadcast  applications  do  not  require  any  gain  uniformity  correction.  The  area  of  HARP  needs  to 
be  increased  significantly  for  x-ray  imaging  applications.  Although  uniform,  large  area  a-Se  films  have  been  developed 
for  direct  FPI,  the  uniformity  requirement  for  HARP  in  indirect  FPI  with  avalanche  gain  is  expected  to  be  much  higher. 
In  this  section  we  will  derive  the  relationship  between  thickness  non-uniformity  and  gain  non-uniformity  for  HARP 
sensors  so  that  the  results  can  be  used  as  a  guideline  for  engineering  of  large  area  HARP  layers  for  FPI. 


During  the  operation  of  HARP,  a  constant  potential  Vb  is  applied  to  the  top  bias  electrode  (ITO).  A  non-uniformity  in  the 
thickness  of  HARP  will  result  in  a  variation  in  ESe  and  hence  a  gain  non-uniformity.  By  substituting  ESe  =  VfJdSe  into  Eq. 
5  and  differentiating  both  sides  of  the  equation,  we  obtain: 

—  A-av  =?&—&-,.  dSe.  (8) 

lngav  gav  dSe  Vh 


Eq.  8  can  be  rewritten  as: 


&av 


^Se  dse 


(9) 


Eq.  9  shows  that  the  percentage  variation  in  avalanche  gain  is  proportional  to  ln(gov)  and  the  percentage  thickness 
variation.  Shown  in  Figure  8  are  the  variations  in  gm,  for  both  positive  and  negative  change  in  dSe  around  8  pm.  Three  Vb 
values  of  840  V,  864  V  and  880  V  were  used,  which  corresponds  to  gav  of  12,  25  and  46,  respectively.  Figure  8  shows 
that  gav  increases  with  decrease  in  dse,  which  is  due  to  an  increase  in  ESe  under  a  constant  Vb.  The  slope  of  the  curve,  i.e. 
percentage  change  in  gav,  is  much  higher  for  negative  change  in  dSe  than  it  is  for  positive.  Hence  when  choosing  the  bias 
potential  Vb  for  a  desired  ESe,  the  smallest  dse  for  the  entire  detector  area  should  be  used,  so  that  the  variation  in  thickness 
is  positive  from  the  reference  thickness.  If  we  assume  the  smallest  dSe  =  8  pm,  the  percentage  changes  in  gav  due  to 
increase  in  dSe  of  1  %  are  19  %,  23  %  and  26  %,  respectively,  for  the  three  Vb  settings  of  840  V,  864  V  and  880  V.  The 
change  in  gav  increases  to  33  %,  40  %,  and  44  %,  respectively,  when  the  variation  in  dSe  increases  to  2  %.  This  range  of 
variation  in  gm,  due  to  thickness  non-uniformity  is  certainly  within  the  capability  of  gain  correction  and  the  dynamic 
range  of  the  electronic  circuits  of  FPI.  With  further  increase  in  dSe  to  4  %,  the  value  for  gav  will  decrease  from  the 
original  setting  of  46,  25  and  12  to  15,  10  and  6,  respectively,  which  is  a  gain  non-uniformity  of  2  to  3  folds.  Although 
the  change  in  gav  is  significant,  it  is  not  expected  to  cause  significant  degradation  in  the  DQE  performance,  as  can  be 
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seen  in  the  example  of  gav  =  46  and  12  shown  in  Figure  6  (a).  Hence  we  cain  conclude  that  a  thickness  non-uniformity  of 
2  %  is  desirable  for  a  HARP  layer,  and  a  non-uniformity  of  4  %  is  still  reasonable. 


Figure  8:  The  avalanche  gain  gm  as  a  function  of  variation  in  the  thickness  of  HARP  under  a  constant  bias  potential  Vb. 

CONCLUSIONS 

The  feasibility  of  a  new  concept  of  indirect  conversion  flat-panel  detector  with  avalanche  gain  is  being  investigated.  The 
avalanche  gain  and  quantum  efficiency  of  an  8  pm  thick  HARP  layer  was  investigated  as  a  function  of  £&.  The 
characteristics  of  HARP  were  applied  to  a  cascaded  linear  system  to  determine  the  operating  conditions  for  x-ray 
imaging.  Our  results  showed  that  by  varying  the  avalanche  gain  between  2  and  46,  x-ray  quantum  noise  limited 
performance  and  a  wide  dynamic  range  can  be  obtained  for  typical  detector  parameters  chosen  for  an  R/F  detector.  Due 
to  the  rapid  increase  in  gm,  as  a  function  of  ESe,  it  is  desirable  to  keep  the  thickness  uniformity  of  a  HARP  layer  within  2 
%,  which  could  keep  the  variation  in  gav  within  ~  40  %. 
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ABSTRACT 


An  indirect  flat-panel  imager  (FPI)  with  avalanche  gain  is  being  investigated  for  low-dose  x-ray 
imaging.  It  is  made  by  optically  coupling  a  structured  x-ray  scintillator  CsI(Tl)  to  an  amorphous 
selenium  (a-Se)  avalanche  photoconductor  called  HARP  (High-gain  Avalanche  Rushing 
Photoconductor).  The  final  electronic  image  is  read  out  using  an  active  matrix  array  of  thin  film 
transistors  (TFT).  We  call  the  proposed  detector  SHARP -AMFPI  (Scintillator  HARP  Active 
Matrix  Flat  Panel  Imager).  The  advantage  of  the  SHARP- AMFPI  is  its  programmable  gain, 
which  can  be  turned  on  during  low  dose  fluoroscopy  to  overcome  electronic  noise,  and  turned  off 
during  high  dose  radiography  to  avoid  pixel  saturation.  The  purpose  of  this  paper  is  to  investigate 
the  important  design  considerations  for  SHARP -AMFPI  such  as  avalanche  gain,  which  depends 
on  both  the  thickness  dSe  and  the  applied  electric  field  ESe  of  the  HARP  layer.  To  determine  the 
optimal  design  parameter  and  operational  conditions  for  HARP,  we  measured  the  Ese 
dependence  of  both  avalanche  gain  and  optical  quantum  efficiency  of  an  8  pm  HARP  layer.  The 
results  were  used  in  a  physical  model  of  HARP  as  well  as  a  linear  cascaded  model  of  the  FPI  to 
determine  the  following  x-ray  imaging  properties  in  both  the  avalanche  and  non-avalanche 
modes  as  a  function  of  (1)  total  gain  (which  is  the  product  of  avalanche  gain  and  optical 
quantum  efficiency);  (2)  linearity;  (3)  dynamic  range;  and  (4)  gain  non-uniformity  resulting  from 
thickness  non-uniformity;  and  (5)  effects  of  direct  x-ray  interaction  in  HARP.  Our  results 
showed  that  a  HARP  layer  thickness  of  8  pm  can  provide  adequate  avalanche  gain  and  sufficient 
dynamic  range  for  x-ray  imaging  applications  to  permit  quantum  limited  operation  over  the 
range  of  exposures  needed  for  radiography  and  fluoroscopy. 

Key  words:  flat-panel  detectors,  avalanche  gain,  amorphous  selenium,  cesium  iodide. 
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1  I.  INTRODUCTION 

2  An  active  matrix  flat-panel  imager  (AMFPI)  based  on  thin  film  transistor  (TFT)  arrays  is  the  most 

3  promising  technology  for  digital  x-ray  imaging  due  to  its  compact  size,  rapid  readout  and  superior 

4  performance  compared  to  screen-films.1  Based  on  the  materials  used  for  x-ray  detection,  AMFPI  are 

5  categorized  as  either  direct  or  indirect  based  on  the  materials  used  for  x-ray  detection  (x-ray 

6  photoconductors  or  scintillators,  respectively).  Direct  AMFPI  have  the  advantages  of  higher  image 

7  resolution  and  a  simpler  TFT  array  structure  that  can  be  manufactured  in  a  standard  facility  for 

8  active  matrix  liquid  crystal  displays  (AMLCD).  Despite  their  differences,  both  types  of  AMFPI  offer 

9  better  image  quality  than  screen-films  and  computed  radiography  (CR). 

10 

1 1  However  before  AMFPI  can  equal  the  performance  of  the  x-ray  image  intensifier  (XRII)  in 

12  fluoroscopy,  two  major  difficulties  must  be  overcome:  (i)  the  ability  to  produce  good  image  quality 

13  at  very  low  dose,  such  as  in  the  dark  part  of  a  fluoroscopy  image  (~0. 1  pR  per  frame),2’3  and  (ii)  the 

14  ability  to  produce  images  at  a  high  frame  rate  without  artifacts,  especially  when  the  radiation 

15  exposure  is  switched  from  radiographic  to  fluoroscopic.4,5  Even  in  well  established  radiographic 

16  applications,  AMFPI  are  facing  new  challenges  because  their  rapid  readout  and  freedom  from 

17  geometric  distortion  have  generated  a  renewed  interest  in  advanced  clinical  applications.  One 

18  example  is  digital  tomosynthesis,6  in  which  a  rapid  sequence  of  images  is  taken  from  different 

19  angles,  and  tomographic  slices  at  different  depths  of  the  patient  anatomy  are  reconstructed.  The 

20  application  of  AMFPI  in  cone-beam  computed  tomography  ’  (CT)  is  also  under  intensive 

21  investigation.  Such  advanced  applications  of  AMFPI  pose  new  challenges  for  the  development  of 

22  the  next  generation  of  detectors,  primarily  regarding  the  ability  of  AMFPI  to  generate  high  quality 

23  images  that  are  quantum  noise  limited  and  free  from  artifacts  at  low  x-ray  exposures  and  high  frame 
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1  rates.  It  has  been  shown  that  the  detective  quantum  efficiency  (DQE)  of  AMFPI  is  inferior  to  XRII 

2  at  the  lower  end  of  the  dose  range  (0.1-1  pR  per  frame)  encountered  in  fluoroscopy.9  In 

3  mammography,  the  detector  is  not  x-ray  quantum  noise  limited  behind  dense  breast,  where  the 

4  average  exposure  is  ~1  mR,  and  the  DQE  of  the  detector  is  severely  compromised  by  the  electronic 

5  noise.  ’  . 

6 

7  The  strategies  for  improving  the  low  dose  performance  of  AMFPI  can  be  divided  into  two  categories: 

8  one  is  to  increase  the  x-ray  to  image  charge  conversion  gain  so  that  the  x-ray  quantum  noise  can 

9  overcome  the  electronic  noise;  and  the  other  is  to  decrease  the  electronic  noise.  These  strategies  are 

10  common  to  both  types  of  AMFPI  since  they  have  approximately  the  same  conversion  gain  and  pixel 

1 1  electronic  noise.1  Several  investigations  are  taking  place  to  increase  the  gain  for  direct  AMFPI  by 

12  developing  other  x-ray  photoconductors.  Promising  results  have  been  shown  for  lead  iodide  (Pbl2) 

13  and  mercuric  iodide  (Hgf).'2  The  conversion  gain  of  these  two  materials  is  5-7  times  higher  than 

14  that  of  current  systems.  The  challenges  for  this  technology  are  to  develop  deposition  techniques  for 

15  large  area  detectors  with  good  imaging  properties,  and  to  ensure  that  both  electrons  and  holes  have 

16  adequate  range  to  reduce  the  probability  of  charge  trapping.13  To  reduce  electronic  noise,  several 

17  investigators  have  proposed  advanced  pixel  designs,  which  incorporate  pixel  amplification  by 

18  adding  at  least  two  more  TFTs  at  each  pixel. 14,15,16  The  electronic  noise  can  thus  be  reduced  to  -700 

19  electrons  per  pixel,  which  is  a  factor  of  2-4  lower  than  existing  AMFPI.  A  challenge  for  this 

20  approach  is  the  practical  implementation  of  the  complex  pixel  design  over  a  large  area  with 

21  consistent  and  uniform  imaging  performance  because  each  pixel  operates  as  an  analog  amplifier 

22  circuit  as  opposed  to  the  simple  switching  device  (digital)  in  existing  AMFPI.  Due  to  these  practical 
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1  challenges  none  of  the  above  approaches  for  increased  gain  and  reduced  electronic  noise  has  been 

2  implemented  in  commercial  large  area  AMFPI. 

3 

4  We  are  proposing  a  new  concept  of  AMFPI  by  adding  avalanche  gain  to  an  indirect  AMFPI,  as 

5  shown  in  Figure  1.  The  optical  photons  emitted  from  the  Csl  layer  upon  absorption  of  x-rays  are 

6  detected  and  amplified  by  an  avalanche  a-Se  photoconductor  called  HARP  (High-gain  Avalanche 

7  Rushing  Photoconductor).  Light  photons  absorbed  by  the  HARP  layer  generate  electron-hole  pairs 

8  near  the  top  interface  (<0.1  pm).  Under  a  sufficiently  high  Ese  (produced  by  applying  positive 

9  voltage  on  the  top  bias  electrode  of  a-Se),  holes  moving  towards  the  bottom  surface  will  undergo 

10  avalanche  multiplication  and  create  more  holes  and  electrons.  The  avalanche  gain  gav  ranges  from 

11  unity  up  to  ~ 1 000  times  depending  on  E$e  and  the  thickness  of  the  a-Se  layer,  dSe-X1  The  electrons 

12  move  to  the  top  bias  electrode  while  holes  are  collected  by  the  pixel  electrodes  and  form  a  charge 

13  image,  which  is  read  out  with  a  two-dimensional  array  of  TFTs,  similar  to  that  in  existing  AMFPI. 

14  The  proposed  detector  is  named  SHARP-AMFPI  (Scintillator-HARP  Active  Matrix  Flat-panel 

1  O  t 

15  Imager).  Rougeot  and  Possin  have  suggested  a  similar  concept,  but  it  has  not  been  implemented. 

16  The  only  other  indirect  AMFPI  approach  involving  an  a-Se  photoconductor  as  an  optical  sensor  did 

1 7  not  have  avalanche  gain. 1 9 

18 

19  Compared  to  existing  indirect  FPI,  SHARP-AMFPI  has  the  following  advantages:  (1)  x-ray  quantum 

20  noise  limited  at  a  much  lower  x-ray  exposure  level;  (2)  better  temporal  performance,  i.e.  lower  lag 

21  and  ghosting  because  of  the  lower  trap  density  in  a-Se  (compared  to  a-Si)  and  the  much  higher 

22  electric  field  2?$e;  (3)  compatibility  with  standard  TFT  array  manufacturing  facility  for  AMLCD;  (4) 

23  ability  to  use  high-resolution  (HR)  type  Csl,  which  cannot  be  used  in  existing  AMFPI  due  to  its 
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1  lower  light  output.  Compared  to  the  high  conversion  gain  x-ray  photoconductors  being  developed 

2  for  direct  FPI,  the  proposed  detector  has  the  following  advantages:  (1)  programmable  gain  by 

3  changing  the  electric  field  Ese  as  opposed  to  the  fixed  high  conversion  gain  of  Hgl2  or  PW2.  One  of 

4  the  practical  problems  of  having  a  fixed  high  gain  is  that  the  large  signal  charge  (at  high  exposure) 

5  requires  a  large  pixel  storage  capacitor  (>15  pF),  which  is  impractical  to  make  especially  for  small 

6  pixel  sizes.  This  means  that  a  detector  that  is  optimized  for  low  dose  (fluoroscopic)  x-ray  imaging 

7  applications  will  not  work  properly  at  high  radiographic  dose  levels.  With  SHARP-AMFPI,  by 

8  appropriately  decreasing  the  ESe  during  high  exposure  applications,  the  signal  charge  is  comparable 

9  to  that  at  low  dose  with  avalanche  gain.  This  simplifies  the  TFT  array  pixel  design  and  permits  a 

10  detector  with  wide  dynamic  range.  (2)  better  temporal  performance  because  both  charge  carriers  in 

1 1  a-Se  have  adequate  range  and  the  probability  of  charge  trapping  is  low.  This  is  in  contrast  to  the  very 

12  short  range  of  one  type  of  carriers  in  Pbl2  (electrons)  and  Hgl2  (holes)  which  results  in  significant 

13  charge  trapping  and  ghosting.  The  advantages  described  above  make  the  proposed  detector  suitable 

14  for  all  x-ray  imaging  modalities  (i.e.  radiography  and  fluoroscopy).  The  purpose  of  this  paper  is  to 

15  investigate  the  important  design  parameters  for  HARP  and  to  establish  the  operational  conditions 

16  that  ensure  adequate  dynamic  range  for  both  fluoroscopy  and  radiography. 

17 

18  Since  avalanche  multiplication  in  a-Se  was  discovered  over  a  decade  ago  the  HARP  video  tube, 

19  which  consist  of  a  HARP  layer  read  out  with  a  scanning  electron  beam,  has  been  commercialized  for 

20  high  sensitivity  and  high  definition  television  applications.  Stable  and  uniform  avalanche 

21  multiplication  has  been  demonstrated  with  very  little  added  noise.20  No  correction  for  fixed  pattern 

22  noise  is  necessary.  The  electric  field  Ese  under  which  avalanche  multiplication  occurs  is  >  70  V/um, 

23  which  is  an  order  of  magnitude  higher  than  that  typically  used  in  direct  AMFPI  incorporating  thick 


4 


1  a-Se  layers.  To  sustain  Ese  of  this  high  magnitude  without  a  significant  increase  in  dark  current  due 

2  to  charge  injection  from  the  bias  electrodes,  blocking  layers  are  required  for  both  electrode 

3  interfaces.  As  shown  in  Figure  1  (b),  the  HARP  sensor  is  a  multi-layer  structure.17  It  contains  a 

4  transparent  indium  tin  oxide  (ITO)  electrode  which  is  biased  positively  during  imaging,  followed  by 

5  a  thin  layer  (~20  nm)  of  cerium  oxide  (CeCf)  and  a  thin  layer  of  lithium  fluoride  (LiF)  doped  a-Se  to 

6  block  injection  of  holes  from  the  ITO.  The  intrinsic  a-Se  photoconductor  is  from  0.5  to  35  pm  thick 

7  depending  on  the  desired  avalanche  gain.  The  bottom  interface  of  the  HARP  structure  is  a  sub- 

8  micron  thin  layer  of  antimony  tri-sulfide  (Sb2S3),  which  is  used  to  block  electron  injection  from  the 

9  negative  bias  electrode  (i.e.  pixel  electrodes  of  the  TFT  array).  Since  the  mobility  of  holes  in  a-Se  is 

10  more  than  30  times  higher  than  that  of  electrons,  only  holes  create  a  significant  number  of  new 

1 1  carriers  by  impact  ionization.  Hence  it  is  essential  for  light  photons  to  enter  from  the  positive  bias 

12  electrode. 

13 

14  Our  previous  calculation  using  a  preliminary  linear  system  model  of  indirect  FPI  with  avalanche 

1 5  gain  predicted  that  an  avalanche  gain  of  gav  =  20  is  needed  to  produce  x-ray  quantum  noise  limited 

16  performance  for  fluoroscopy  and  low  dose  mammography  (tomosynthesis).21  This  requires  a  HARP 

17  layer  thickness  of  4  pm  or  above.  For  the  present  paper  we  investigated  the  effects  of  the  optical 

1 8  imaging  properties  of  HARP  on  the  x-ray  imaging  performance  of  SHARP- AMFPI,  so  that 

19  appropriate  design  parameters  and  operational  conditions  could  be  determined.  We  measured  the  Ese 

20  dependence  of  both  avalanche  gain  and  optical  quantum  efficiency  of  an  8  pm  HARP  layer.  The 

21  results  were  used  in  a  physical  model  of  HARP  as  well  as  a  linear  cascaded  model  of  the  FPI  to 

22  determine  the  following  x-ray  imaging  properties  in  both  the  avalanche  and  non-avalanche  modes  as 

23  a  function  of  dse  and  £&:  (1)  total  gain  (which  is  the  product  of  avalanche  gain  and  optical  quantum 
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1  efficiency);  (2)  linearity;  (3)  dynamic  range;  (4)  gain  non-uniformity  resulting  from  thickness  non- 

2  uniformity;  and  (5)  the  effects  of  direct  x-ray  interaction  in  HARP. 

3 

4  II.  PHOTOSENSITIVITY  OF  HARP 

5  The  x-ray  imaging  performance  of  indirect  FPI  with  avalanche  gain  relies  on  the  photosensitivity  of 

6  HARP.  Understanding  the  dependence  of  photosensitivity  on  Ese  and  the  layer  thickness  dSe  is 

7  crucial  for  guiding  the  choice  of  HARP  design  parameters  and  operational  conditions  for  both  the 

8  avalanche  and  non-avalanche  modes.  In  this  section,  we  will  first  present  the  experimental  method 

9  and  results  of  the  HARP  photosensitivity  measurements,  and  then  describe  the  theory  of  £& 

10  dependence  of  photosensitivity,  which  is  proportional  to  both  the  optical  quantum  efficiency  (QE)  r\ 

1 1  and  the  avalanche  gain  gav.  The  theory  and  experimental  data  will  then  be  used  to  determine  r|  and 

12  gav  as  a  function  of  E$e  in  the  context  of  x-ray  imaging. 

13 

14  A.  Experimental  Measurement  of  Photosensitivity 

15  1.  Method 

16  The  experimental  apparatus  for  photosensitivity  measurement  is  shown  in  Figure  2.  A  grating 

17  monochromator  was  used  to  generate  optical  photons  of  a  single  wavelength  X.  The  values  of  X  used 

18  in  our  experiments  were  in  the  ranged  of  380-600  nm.  The  intensity  of  the  output  beam  was 

19  attenuated  by  neutral-density  (ND)  filters.  A  beamsplitter  was  used  to  direct  a  fraction  (90  %)  of  the 

20  attenuated  beam  to  a  silicon  (Si)  photodiode  (PD),  which  was  used  to  monitor  the  beam  intensity 

21  during  the  experiment.  The  optical  photons  passing  directly  through  the  beam  splitter  (10  %)  were 

22  detected  by  the  HARP  video  tube.  Before  photosensitivity  measurements  the  input  light  intensity  to 
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1  the  HARP  was  calibrated  such  that  a  constant  input  power  was  used  for  each  X.  This  was  achieved 

2  by  adjustment  of  the  ND  filters  based  on  the  light  intensity  monitored  by  the  Si  PD. 

3 


4  For  each  X  and  Ese,  the  signal  current  S  from  the  HARP  tube  was  recorded.  The  photosensitivity  was 

5  quantified  as  the  number  of  signal  charge  carriers  generated  per  absorbed  light  photon.  This  is 

6  referred  to  as  the  effective  quantum  efficiency  T]*,  and  it  is  the  product  of  r\  and  gav.  The  values  of  r|* 

7  were  experimentally  determined  from  the  measurement  of  HARP  tube  current  S  using: 
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«  S/e 
~  IT/hv 


(2) 


9  where  e  is  the  electron  charge,  /  is  the  HARP  input  light  intensity  power  given  in  Watts  (W),  h  vis 

10  the  energy  of  each  incident  photon,  and  T  is  the  fraction  of  the  photons  incident  on  the  HARP  tube 

1 1  reaching  the  a-Se  layer.  The  correction  factor  T  represents  the  light  attenuation  by  the  glass  faceplate, 

12  the  ITO  bias  electrode  and  the  CeC>2  blocking  layer  for  HARP.  The  value  of  T  is  a  function  of  the 

1 3  wavelength  X  and  was  experimentally  measured. 

14 

15  2.  Measurement  results: 

16  The  measured  effective  quantum  efficiency  r\*  is  shown  in  Figure  3  as  a  function  of  £&  for  an  8  pm 

17  thick  HARP  layer.  For  blue  light  (A=400  nm),  r\*  reaches  a  plateau  at  ESe  >  20  V/pm.  This 

1 8  corresponds  to  a  saturation  of  optical  quantum  efficiency  r|  at  unity.  With  further  increase  in  Ese  to  > 

19  80  V/pm,  rf  starts  to  increase  again,  which  shows  the  onset  of  avalanche.  However  for  green  light, 

20  the  Ese  dependence  of  q*  is  quite  different.  As  shown  in  Figure  3  for  A=540  nm,  which  is  the  peak  of 

21  the  emission  spectrum  of  Csl  (Tl),  r\*  starts  at  a  much  smaller  value  compared  to  that  for  blue  light. 

22  It  then  increases  continuously  with  Ese  without  reaching  a  plateau  before  avalanche  begins  at  ~80 
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1  V/|im.  This  is  because  the  optical  q  for  green  light  is  lower  and  does  not  reach  unity  until  Ese  is  well 

2  into  the  range  of  avalanche  multiplication.  To  choose  the  proper  operational  condition  for  HARP  in 

3  different  x-ray  imaging  applications,  we  need  to  understand  separately  the  Ese  dependence  of  r|  and 

4  gav.  This  will  be  discussed  in  section  1IB. 

5 

6  B.  Determination  of  Quantum  Efficiency  and  Avalanche  Gain 

7  The  effective  quantum  efficiency  q*  of  HARP  is  proportional  to  both  q  and  gav.  The  value  of  q  in  a- 

8  Se  depends  on  both  ESe  and  the  wavelength  X  of  the  incident  photon.  It  was  found  by  Pai  and  Enck 

9  that  the  Ese  and  X  dependence  of  photogeneration  in  a-Se  follow  the  Onsager  mechanism.22  Onsager 

10  theory  states  that  every  absorbed  photon  creates  a  pair  of  thermalized  carriers  bound  by  their 

1 1  Coulombic  attraction.  The  initial  separation  ro  of  the  charge  pair  increases  with  incident  photon 

12  energy.  The  probability  of  dissociation  for  the  bound  pairs  into  free  carriers  increases  with  r0  and  the 

13  applied  external  field  E&  with  the  relationship: 

14  eESer0  kT  />.’ 

e2  1 

with  A  = - 

4/T££0kT 

15  where  4»0  is  defined  as  the  efficiency  of  production  of  thermalized  charge  pairs  per  absorbed  photon 

16  and  is  independent  of  £&.  Pai  and  Enck  showed  that  the  room  temperature  measurement  of  q  can  be 

17  best  fitted  with  the  expression  in  Eq.  1  using  4>0=1  and  ro  in  the  range  between  0.8  -  7.0  nm  for  X 

18  values  between  620  and  400  nm.  For  the  rest  of  the  discussion  in  this  paper,  we  will  focus  on  A=540 

19  nm,  which  is  the  peak  of  the  emission  spectrum  for  Csl  (Tl).  Shown  in  Figure  4  are  the  values  of  q 

20  calculated  using  Eq.  1  and  ro  =  1.7  nm.  Also  plotted  in  Figure  4  are  the  q*  measured  at  the  same 
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1  wavelength  for  dse  =  8  pm.  These  results  show  that  the  Onsager  theory  provides  a  reasonable 

2  prediction  for  the  Ese  dependence  of  r\*  in  HARP  for  Ege  <  70  where  gav=  1  •  The  calculated  r|  using 

3  Onsager  theory  was  then  divided  into  the  measured  rf  so  that  the  Ese  dependence  of  avalanche  gain 

4  gav  can  be  obtained.  The  results  are  shown  in  Figure  5  (a),  where  gav  =  46  at  Ese  =  1 10  V/pm.  Figure 

5  5  (a)  shows  that  gav  increases  very  rapidly  with  £&>,  hence  choosing  the  appropriate  E$e  for  HARP  is 

6  crucial.  To  determine  the  x-ray  response  of  SHARP -AMFPI  it  is  necessary  to  understand  the 

7  continuous  dependence  of  gav  on  Ese  because  as  image  charge  accumulates  on  the  pixels,  the  voltage 

8  drop  across  the  HARP  layer  will  decrease.  This  change  in  ESe  will  produce  a  corresponding  change 

9  in  gav  and  affect  the  x-ray  sensitivity,  linearity  and  dynamic  range  of  the  FPL 
10 

1 1  The  avalanche  gain  characteristics  in  HARP  have  been  investigated  as  a  function  of  dse  and  ESe23 

12  According  to  the  theory  of  avalanche  multiplication,24  when  Ese  across  the  thickness  dse  is  uniform, 

13  and  the  distance  through  which  avalanche  occurs  is  equal  to  dse  (which  is  true  for  an  8  pm  HARP 

14  because  photons  are  absorbed  within  0. 1  pm  at  the  entrance  side  of  the  a-Se  layer),  gav  is  given  by: 

Sav  p-cced*(p-a)  V) 

16  where  a  and  P  are  the  impact  ionization  coefficients  of  electrons  and  holes,  respectively.  The  value 

17  of  p  is  much  greater  than  a  since  in  a-Se  the  drift  mobility  of  holes  (0. 1 8  cm2/V s)  is  much  higher 

18  than  that  of  electrons  (0.003  cm2/Vs).  This  is  why  in  HARP  light  always  enters  from  the  positive 

19  bias  electrode  to  maximize  avalanche  gain.  Under  the  assumption  that  the  contribution  of  electrons 

20  to  avalanche  gain  is  negligible,  Eq.  2  can  be  simplified  to: 

21  g„=e^e  (3) 
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1  with  p  given  by: 


-Pi 

2  fi  =  /3xeE*.  (4) 

3  The  values  for  Pi  and  P2  can  be  determined  by  curve  fitting  of  experimental  measurement  of  gav.  By 

4  taking  the  logarithms  of  both  sides  of  Eq.  3  twice,  we  obtain: 

5  HlngJ  =  inWs,)~.  (5) 

Ese 

6  which  indicates  that  the  quantity  ln(lngav)  is  linearly  related  to  1  /Ese.  Shown  in  Figure  5  (b)  is  the 

7  experimentally  determined  gav  data  from  Figure  5(a)  replotted  in  the  form  of  ln(lngav)  as  a  function 

8  of  1  JEse-  The  best  linear  fit  to  the  experimental  data  is  also  shown  in  Figure  5  (b).  From  the  intercept 

9  and  the  slope  of  the  fitted  line  we  obtained  Pi=5.5  x  103  pm'1  and  p2  =  1.029  x  103  V/pm.  Using 

10  these  pi  and  P2  values  and  the  analytical  expression  for  gav  shown  in  Eqs.  3  and  4,  the  fitted  gav 

1 1  values  were  calculated  and  are  shown  in  Figure  5  (a).  It  shows  that  the  calculated  gav  using  the  fitted 

1 2  parameters  agrees  well  with  experimentally  determined  gav  for  ESe  <110  V/pm.  Hence  the  fitted 

13  parameters  Pi  and  p2  will  be  used  in  the  following  investigation  of  the  potential  x-ray  imaging 

14  performance  of  SHARP- AMFPI. 

15 

1 6  III.  POTENTIAL  X-RAY  IMAGING  PERFORMANCE 

17  A.  Cascaded  Linear  System  Model  for  SHARP-AMFPI 

1 8  Previously  we  have  developed  a  cascaded  linear  system  model  for  indirect  FPI  with  avalanche 

19  gain.21  This  preliminary  investigation  did  not  include  K- fluorescence  reabsorption  in  Csl  or  direct  x- 

20  ray  interaction  in  HARP.  We  showed  that  the  DQE(/)  is  related  to  the  avalanche  gain  gav  through: 


10 


DQE(/)  =  ■ 


^  T1  2  /  At2/  S\  ,  (2g 

it1  2 


wwh  -  &av  - 


Tb\f)Tc\f)Ta\f) 
f  Vs°A)]T2{f)+ 


(<*P  gavgth<xn)  % 


2  where  (70  is  the  number  of  incident  x-ray  photons  per  unit  area,  Tb(f),  TC(J),  and  Ta(f)  are  the  MTFs 

3  due  to  blur  in  Csl,  the  optical  coupling  between  Csl  and  HARP,  and  the  pixel  aperture  function, 

4  respectively.  As,  v\x,  and  g  are  the  Swank  factor,  x-ray  quantum  absorption  efficiency  and  x-ray  to 

5  photon  conversion  gain  of  the  Csl  layer,  respectively.  The  factors  related  to  the  TFT  array  are  the 

6  electronic  noise  S„  and  the  pixel  width  aP.  The  denominator  of  Eq.  6  is  the  dose  normalized  NPS. 

7  The  first  term  is  the  x-ray  quantum  noise,  the  second  term  is  the  secondary  quantum  noise  associated 

8  with  the  variance  in  gain  in  the  conversion  and  the  avalanche  process,  and  the  third  term  is  the 

9  contribution  from  the  electronic  noise  of  the  TFT  array.  This  expression  shows  that  the  electronic 

10  noise  is  more  important  at  low  dose,  and  is  inversely  proportional  to  gav2.  Therefore  introducing 

1 1  avalanche  gain  is  an  effective  method  for  reducing  noise  at  low  dose.  The  secondary  quantum  noise 

1 2  for  gav »  1  is  approximately  twice  that  for  gav=  1,  assuming  r\  is  the  same  in  both  cases.  This 

13  excess  secondary  quantum  noise  is  the  penalty  for  avalanche  gain.  However  for  SHARP-AMFPI 

14  with  high  conversion  gain  g  and  efficient  optical  coupling  between  Csl  and  HARP,  i.e.  grixarj  »  2 , 


15  the  secondary  quantum  noise  is  negligible  compared  to  the  x-ray  quantum  noise.  Hence  the  DQE  of 

16  SHARP-AMFPI  will  not  be  degraded  by  the  noise  associated  with  avalanche  gain.  This  has  been 

17  demonstrated  experimentally  by  coupling  a  HARP  tube  to  the  output  of  an  x-ray  image  intensifier 

1 8  (XRII)  and  measuring  the  DQE  as  a  function  of  avalanche  gain  gav.25  It  was  shown  that  DQE 


19  remains  essentially  constant  as  a  function  of  gav  (with  electronic  noise  negligible  in  all  cases). 

20 
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1  In  a  later  investigation  of  the  inherent  imaging  performance  of  Csl  layers,  we  included  the  effect  of 

2  K-fluorescence  reabsorption  of  Csl  in  the  cascaded  linear  system  model  for  indirect  AMFPI.26  Our 

3  results  showed  that  the  K-fluorescence  reabsorption  in  Csl  has  two  effects:  (1)  degradation  of  the 

4  Swank  factor  for  x-ray  energies  above  the  K-edges;  and  (2)  additional  image  blur,  which  accounts 

5  for  a  drop  in  presampling  MTF  of  20  -  30  %  at  10  cycles/mm  for  Csl  thickness  between  150  -  600 

6  pm.  For  the  present  paper,  a  more  complete  cascaded  linear  system  model  was  developed  that 

7  includes  the  effect  of  direct  x-ray  interaction  in  HARP.  Shown  in  Figure  6  are  the  serial  and  parallel 

8  cascaded  stages  involved  in  the  formation  of  an  x-ray  image  by  SHARP-AMFPI.  Each  stage  can  be 

9  categorized  as  one  of  the  following  six  processes:  (1)  gain  or  selection,  (2)  stochastic  blurring,  (3) 

10  deterministic  blurring,  (4)  aliasing,  (5)  addition,  and  (6)  parallel  processes.  There  are  two  parallel 

1 1  processes:  one  associated  with  the  probability  of  K-fluorescence  generation  and  reabsorption  in  Csl, 

12  and  the  other  with  direct  interaction  of  x-rays  in  the  HARP  layer.  We  have  previously  described  in 

13  detail  the  transfer  of  spatial  frequency  dependent  signal,  ® (/),  and  NPS,  S(f)  through  all  the  stages 

14  except  direct  x-ray  interaction  in  the  HARP  layer.  This  process  will  be  described  in  Section  C  of  this 

15  paper. 

16 

17  Using  the  cascaded  linear  system  model  shown  in  Figure  6,  we  investigated  two  important  aspects  of 

18  the  x-ray  imaging  performance  of  SHARP-AMFPI.  One  is  the  pixel  x-ray  response  of  the  detector, 

19  i.e.  the  linearity  and  dynamic  range,  for  both  the  avalanche  and  non-avalanche  operational  modes; 

20  the  other  is  the  effect  of  direct  x-ray  interaction  in  HARP  on  the  DQE  of  the  detector.  These  will  be 

21  described  in  Sections  B  and  C,  respectively. 

22 

23  B.  Pixel  X-ray  Response 
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In  this  section,  we  will  apply  the  £&>  dependence  of  gav  and  r)  derived  previously  in  Section  II  (B)  to 
the  cascaded  linear  system  model  to  determine  the  x-ray  response  of  SHARP-AMFPI.  Due  to  the 
difference  in  resolution  and  avalanche  gain  requirement  for  mammography  and 
radiography/fluoroscopy  (R/F),  the  calculation  will  be  performed  separately  for  these  two  imaging 
applications.  From  the  cascaded  linear  system  model  the  mean  signal  from  each  pixel  of  the  detector 
is  given  by: 


0  _  aP2angqVEab 
s  W 


*  (7) 


where  W  is  the  energy  required  to  generate  an  optical  photon  in  Csl,  and  Eab  is  the  x-ray  energy 
absorbed  by  Csl  per  unit  area.  E„b  includes  contributions  from  all  three  parallel  paths  associated  with 
K-fluorescence: 

Eab^q0(E)rix(E)[(l-Pk)E  +  (E~Ek)Pk+EkPkfk]  (8) 

E 

where  qo(E)  is  the  incident  x-ray  photons  per  unit  area  as  a  function  of  energy  E  for  a  given  x-ray 
spectrum,  Pkis  the  probability  of  K-fluorescence,  Ek  is  the  energy  of  the  K-fluorescence  photon,  and 
fk  is  the  fraction  of  K-fluorescence  reabsorbed  in  the  Csl  layer.  Each  of  the  K-fluorescence  related 
terms  in  Eq.  8  were  determined  separately  for  the  Ka  and  Kp  photons  from  both  Cs  and  I  atoms  and 
integrated  to  obtain  the  final  result  of  Eab. 


Our  previous  investigation  of  the  inherent  x-ray  imaging  performance  of  Csl  (Tl)  has  shown  that  the 
high  resolution  (HR)  type  has  essentially  no  depth  dependent  blur  and  as  a  result  provides  better 
DQE(/)  at  high  spatial  frequencies  compared  to  the  Csl  optimized  for  high  light  output  (HL).26  This 
was  achieved  at  the  cost  of  a  lower  light  output  (60  %)  for  the  HR  compared  to  the  HL  Csl.  For  the 
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1  calculation  in  this  paper,  we  used  150  pm  HR  Csl  for  mammography  and  600  pm  HL  Csl  for  R/F 

2  applications.  Our  previously  measured  values  for^s,  W  and  presampling  MTF  of  these  Csl  layers 

3  were  used  in  the  cascaded  linear  system  model.  Since  we  have  not  yet  implemented  the  depth 

4  dependent  blur  of  Csl  in  our  cascaded  linear  system  model,  the  DQE(/)  at  high  frequencies  for  the 

5  600  pm  HL  Csl  layer  (for  R/F)  will  be  overestimated.  Table  I  summarizes  the  detector  parameters 

6  and  HARP  operating  conditions  chosen  for  each  application.  For  mammography,  we  consider  two 

7  different  image  acquisition  modes:  screening  mammography  and  tomosynthesis  (where  several 

8  images  of  the  breast  are  acquired  from  different  viewing  angles  for  tomographic  reconstruction  of 

9  3D  images). 

10 

1 1  1.  Radiography/Fluoroscopy 

12  Shown  in  Figure  7  (a)  is  the  DQE(/)  calculated  using  the  cascaded  linear  system  model  for  the 

13  detector  with  different  gav  settings.  The  calculation  was  performed  for  the  lowest  exposure  in 

14  fluoroscopy  (0. 1  pR),  where  the  detector  performance  is  most  susceptible  to  electronic  noise.  It 

15  shows  that  gav  =  46  is  sufficient  to  bring  the  DQE(/)  to  its  theoretical  limit  (i.e.  a  detector  without 

16  electronic  noise).  The  DQE(/)  for  the  next  highest  gav  =  12  measured  at  ESe  =  105  V/pm  also 

17  provides  reasonable  DQE.  Thus  an  ESe  between  105  and  1 10  V/pm  will  be  adequate  for  fluoroscopy. 

1 8  Figure  7  (b)  shows  the  DQE(/)  calculated  for  the  lowest  radiographic  exposure  of  30  pR.  It  shows 

19  that  ESe  =  95  V/pm  can  provide  sufficient  q  (0.3 1)  and  gav  (2.3).  Hence  the  Ese  chosen  in  Table  I  is 

20  adequate  for  the  dual-mode  operation  of  a  R/F  detector. 

21 

22  Shown  in  Figure  8  are  the  calculated  pixel  x-ray  responses  for  the  R/F  detector  under  two  different 

23  operational  conditions.  Here  we  assumed  that  the  dominant  pixel  capacitance  is  the  capacitance  of 
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1  the  HARP  layer.  This  condition  results  in  the  largest  drop  in  the  potential  across  the  HARP  layer  for 

2  a  given  x-ray  exposure,  and  creates  the  biggest  effect  on  linearity.  Figure  8  (a)  shows  that  for 

3  fluoroscopy,  the  initial  £'*=110  V/pm  provides  a  gav  of  46.  Within  the  regular  exposure  range  of 

4  0. 1  -  1 0  pR  per  frame  in  fluoroscopy,  the  detector  has  good  linearity,  i.e.  an  essentially  constant  gav. 

5  With  further  increase  in  exposure,  which  corresponds  to  the  exposure  near  or  beyond  the  edge  of  the 

6  human  body,  the  x-ray  response  becomes  sublinear  with  a  decrease  in  gav ■  This  decrease  in  gav 

7  ensures  a  wide  dynamic  range  for  the  detector  because  there  is  no  “hard”  saturation.  The  pixel 

8  potential  will  continue  to  increase  as  a  function  of  x-ray  exposure  until  it  reaches  the  predetermined 

9  Vmax  for  high  voltage  protection  of  the  TFTs.  The  method  for  high  voltage  protection  was  described 

10  in  detail  elsewhere.27.  The  pixel  potential  for  the  exposure  of  1  mR  in  the  fluoroscopy  mode  is  Vp  = 

11  56  V,  which  is  safe  for  a-Si  TFTs.  Figure  8  (b)  shows  that  with  a  much  lower  gav  =  2.3  programmed 

12  for  radiography,  the  detector  response  for  the  regular  exposure  range  of  30  pR  -  3  mR  is  essentially 

13  linear  with  very  little  change  in  gav.  Comparing  Figures  8  (a)  and  (b),  the  image  charge  on  the  pixel 

14  is  in  the  same  range  despite  their  large  difference  in  exposure  range.  This  demonstrates  the  wide 

15  dynamic  range  (over  5  orders  of  magnitude)  of  the  detector  that  is  made  possible  by  varying  gav  for 

16  different  x-ray  imaging  applications.  This  is  an  advantage  over  detectors  with  constant  high  gain, 

17  which  have  the  problem  of  pixel  saturation  for  radiographic  applications. 

18 

19  2.  Mammography 

20  Figure  9  (a)  and  (b)  shows  the  DQE  calculated  for  mammography  using  the  detector  parameters  and 

21  operating  conditions  listed  in  Table  I  for  tomosynthesis  and  screening  applications.  Figure  9  shows 

22  that  the  choice  of  ESe  for  both  operational  modes  is  adequate  to  achieve  the  theoretical  limit  of  DQE 

23  without  electronic  noise. 
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2  The  pixel  x-ray  response  for  mammography  was  calculated  using  Eq.  7  and  the  results  are  shown  in 

3  Figures  10  (a)  and  (b)  for  tomosynthesis  and  screening  mammography,  respectively.  Two  sets  of 

4  curves  are  shown  in  each  graph  of  Figure  10,  one  without  additional  pixel  storage  capacitance,  and 

5  the  other  with  a  storage  capacitor  of  Cst  -  0.5  pF.  The  minimum  exposure  in  each  graph  of  Figure  10 

6  corresponds  to  the  detector  exposure  behind  dense  breast  tissue,  and  the  maximum  corresponds  to 

7  raw  exposure.  Figure  10  shows  that  when  there  is  no  additional  storage  capacitance,  the  signal 

8  increases  sublinearly  as  a  function  of  the  x-ray  exposure.  This  is  because  the  pixel  capacitance  due 

9  to  the  F1ARP  layer  for  mammography  is  only  -0.03  pF  with  ap  =  70  pm,  which  is  much  smaller  than 

10  that  for  the  R/F  detector  (0.28  pF).  Hence  the  signal  accumulated  on  each  pixel  creates  a  much  more 

1 1  significant  drop  in  ESe,  which  causes  a  rapid  decrease  in  gav.  The  maximum  image  charge  shown  in 

12  Figures  10  (a)  and  10  (b)  corresponds  to  a  pixel  potential  of -100  and  210  V,  respectively.  These 

13  potentials  exceed  the  maximum  safe  voltage  for  the  TFTs.  One  method  for  alleviating  this  problem 

14  while  maintaining  the  same  dynamic  range  is  to  add  a  pixel  storage  capacitance  Cst.  With  Cst  =  0.5 

15  pF,  the  pixel  potential  rises  much  more  slowly  and  the  x-ray  response  of  the  detector  becomes  more 

16  linear,  as  shown  in  Figure  10.  The  pixel  potential  corresponding  to  the  maximum  exposures  in 

17  Figures  10  (a)  and  (b)  is  only  25  and  50  V,  respectively,  both  of  which  are  safe  for  the  TFTs. 

18 

19  C.  Direct  X-ray  Interaction  in  HARP 

20  In  indirect  AMFPI  using  phosphor  and  a-Si  photodiodes,  direct  interaction  of  x-rays  in  the 

21  photodiodes  produces  additional  noise.28  This  is  because  the  x-ray  to  charge  conversion  gain  in  the 

22  photodiodes  is  45  times  higher  than  that  in  the  phosphors,  and  therefore  each  direct  interaction  x-ray 

23  produces  much  higher  signal  (and  noise)  than  that  absorbed  in  the  phosphor.  For  our  proposed 
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1  detector,  the  HARP  sensor  (e.g.  dse  =  8  pm)  is  much  thicker  than  the  a-Si  photodiodes  (~  1  pm)  in 

2  existing  indirect  FPI,  hence  there  will  be  more  direct  x-ray  interaction  events.  Fortunately  the  x-ray 

3  to  charge  conversion  gain  in  HARP  (Wse  -  16  eV  at  ESe  =  100  V)  is  very  similar  to  that  in  Csl,  thus 

4  direct  x-ray  interaction  in  HARP  is  not  expected  to  add  additional  noise  at  low  spatial  frequencies. 

5  At  high  spatial  frequencies,  however,  the  effect  of  direct  x-ray  interaction  could  be  more  significant. 

6  This  is  because  both  the  signal  and  the  NPS  due  to  x-rays  absorbed  in  Csl  decrease  as  a  function  of 

7  spatial  frequency  due  to  blur  in  Csl,  whereas  the  signal  and  NPS  due  to  direct  x-ray  interaction  in 

8  HARP  are  essentially  white  (without  including  the  effect  of  K-fluorescence).  Furthermore,  when 

9  HARP  is  operated  in  the  avalanche  mode,  the  x-rays  absorbed  by  HARP  will  experience  different 

10  avalanche  gain  depending  on  their  depth  of  interaction.  As  shown  in  Figure  11,  x-rays  absorbed  in 

1 1  HARP  at  a  depth  x  from  the  top  (light  entrance)  interface  have  an  effective  thickness  for  avalanche 

12  multiplication  of  ( dse  -  x).  This  could  lead  to  a  depth  dependent  variation  in  gav,  and  add  additional 

13  noise.30  In  the  following  we  will  derive  the  propagation  of  signal  and  noise  due  to  direct  x-ray 

14  interaction  in  HARP  using  the  flow  diagram  shown  in  Figure  6,  and  investigate  its  effect  on  the 

1 5  DQE  of  SHARP-AMFPI. 

16 

17  1.  Signal  and  noise  propagation  for  direct  x-ray  interaction 

18  As  shown  in  Figure  6,  the  input  to  the  parallel  path  for  direct  x-ray  interaction  in  HARP  is  the  x-rays 

19  transmitted  through  the  Csl  layer,  qi,  which  is  a  function  of  x-ray  energy.  For  clarity  we  assume  a 

20  monoenergetic  x-ray  beam  in  the  following  derivation  of  theory,  although  we  used  the  full  x-ray 

21  spectrum  in  our  calculation  of  direct  x-ray  interaction.  We  ignored  the  effect  of  K-fluorescence 

22  reabsorption  in  HARP.  Under  this  assumption,  both  the  signal  and  NPS  for  direct  x-ray  interaction 

23  are  white  prior  to  the  pixel  aperture  function,  hence  we  could  omit  the  frequency  dependence  in  the 
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following  derivation.  The  x-ray  signal  produced  at  depth  x  into  the  HARP  layer  by  a  small  thickness 
dx  of  a-Se  is  given  by: 


<P(x)  =  qxe  MXgav(x)gSeH dx »  (9) 

where  p  is  the  linear  attenuation  coefficient  of  a-Se,  and  gse  is  the  average  x-ray  to  charge 
conversion  gain  which  is  equal  to  E/Wse.  The  depth  dependence  in  gav  can  be  derived  from  Eq.  3  as 
gav(x)  =  exp [j3(dSe  -x)] .  The  NPS  due  to  direct  interaction  as  a  function  of  depth,  S(x),  is  given 


by:30 

S(x)  =  q{  e~ux[(gSe2  +  a\sr  )gj(x)  +  gSe°2gm  (x)]judx .  (10) 

Integrating  both  Eqs.  9  and  10  with  respect  to  x  from  0  to  dse,  we  obtain  the  total  signal  and  NPS  due 
to  direct  x-ray  interaction: 


^  ...  (hMgsei^3' 

M  +  P 

and: 


(ID 


S  = 


<lMgse2+<r2gJ(e 


2/M* 


H  +  '1/3 


(12) 


2.  Calculation  results 

Using  Eqs.  1 1  and  12,  we  calculated  the  contribution  to  the  total  signal  and  NPS  by  the  direct 
interaction  of  x-rays  in  HARP  using  the  parameters  shown  in  Table  I  for  an  R/F  detector.  With  a 
RQA5  x-ray  spectrum  and  a  600  pm  HL  layer  of  Csl,  the  fraction  of  x-rays  transmitted  through  the 
Csl  layer  is  13.2%,  out  of  which  ~1  %  is  absorbed  by  the  HARP  layer.  Shown  in  Figure  12  (a)  is  a 
comparison  between  the  presampling  signal  spectra  (before  the  pixel  aperture  function)  due  to  x-rays 
absorbed  in  Csl  and  direct  x-ray  interaction  in  HARP.  At  low  frequencies,  the  signal  due  to  direct  x- 
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1  ray  interaction  is  less  than  0.2  %  of  that  due  to  x-ray  absorbed  in  Csl.  This  is  consistent  with  the 

2  difference  in  x-ray  quantum  efficiency  between  Csl  and  HARP.  At  high  frequencies  (e.g.  5 

3  cycles/mm),  the  signal  from  Csl  decreases  due  to  the  blur  in  Csl,  however  the  signal  from  direct  x- 

4  ray  interaction  in  HARP  stays  constant.  But  the  latter  is  still  negligible.  Shown  in  Figure  12  (b)  is  the 

5  NPS  (before  pixel  aperture  function)  comparison.  At  5  cycles/mm,  the  NPS  due  to  direct  x-ray 

6  interaction  becomes  significant  (~17  %  of  the  NPS  due  to  Csl).  However  below  the  Nyquist 

7  frequency  of  2.5  cycles/mm,  it  is  still  negligible  (<  3  %).  This  is  why  there  is  no  significant  change 

8  in  DQE  after  including  the  direct  x-ray  interaction  events  in  HARP,  as  shown  in  Figure  13. 

9  Therefore  we  can  conclude  that  direct  x-ray  interaction  in  HARP  in  an  indirect  FPI  with  avalanche 

10  gain  has  no  degradation  effect  on  DQE  of  the  detector. 

11 

1 2  IV.  EFFECT  OF  HARP  THICKNESS  UNIFORMITY 

13  An  important  issue  for  making  a  large  area  HARP  is  its  thickness  uniformity.  Currently  HARP  films 

14  can  be  made  very  uniform  for  an  area  of  1-2”  in  diameter  with  essentially  no  visible  gain  variation  in 

15  the  images.  The  images  currently  produced  by  HARP  for  broadcast  applications  do  not  require  any 

16  gain  uniformity  correction.  The  area  of  HARP  needs  to  be  increased  significantly  for  x-ray  imaging 

17  applications.  Uniform,  large  area  a-Se  films  have  already  been  developed  for  direct  FPI,  however  the 

18  uniformity  requirement  for  HARP  in  indirect  FPI  with  avalanche  gain  is  expected  to  be  much  higher. 

19  In  this  section  we  will  derive  the  relationship  between  thickness  non-uniformity  and  gain  non- 

20  uniformity  for  HARP  sensors  so  that  the  results  can  be  used  as  a  guideline  for  engineering  of  large 

21  area  HARP  layers  for  FPI. 

22 

23  A.  Relationship  between  HARP  Avalanche  Gain  and  Thickness  Uniformity 
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1  During  the  operation  of  HARP,  a  constant  potential  Vb  is  applied  to  the  top  ITO  bias  electrode,  as 

2  shown  in  Figure  1  (b).  After  each  image  readout,  the  TFT  array  resets  the  potential  on  pixel 

3  electrodes  to  Vp  =  0  V.  A  non-uniformity  in  the  thickness  of  HARP  will  result  in  a  variation  in  Ese 

4  and  hence  a  spatial  non-uniformity  in  gain.  An  increase  in  dse  under  a  constant  Vb  causes  a  decrease 

5  in  Ese ■  Although  on  one  hand  an  increase  in  dse  tends  to  increase  gav  (Eq.  3),  it  also  results  in  a 

6  simultaneous  decrease  in  Ese  under  a  constant  F>.  The  overall  effect  is  a  decrease  in  gav  because  of 

7  its  additional  exponential  dependence  on  ESe  shown  in  Eq.  4.  By  substituting  ESe  -  Vb/dSe  into  Eq.  5 

8  and  differentiating  both  sides,  we  obtain: 
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-  —  ^  ^  -  Aa  d. 

^gm  gav  dSe  Vb  Se 


(13) 


10  Eq.  13  can  be  rewritten  as: 

11 


(14) 


12  Eq.  14  shows  that  with  p2  >  -E&,  which  is  true  for  HARP,  a  positive  change  in  dse  results  in  a 

13  negative  change  in  gav.  Eq.  14  also  shows  that  the  percentage  variation  in  avalanche  gain  is 

14  proportional  to  ln(gav)  and  the  percentage  thickness  variation,  hence  HARP  layers  operated  with 

1 5  higher  gav  are  more  susceptible  to  gain  non-uniformity  caused  by  thickness  non-uniformity. 

16 

17  B.  Calculation  of  Gain  Non-uniformity 

18  Using  Eqs.  3  and  4,  and  the  values  for  pi  and  p2  determined  in  Section  II,  the  variations  in  gav  were 

1 9  calculated  for  both  positive  and  negative  change  in  dse  around  8  pm,  and  results  are  shown  in  Figure 

20  14.  Three  Vb  values  of  840  V,  864  V  and  880  V  were  used,  which  corresponds  to  gav  of  12,  25  and 

21  46,  respectively.  Figure  14  shows  that  the  slope  of  the  curve,  i.e.  percentage  change  in  gav,  is  much 
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1  higher  for  negative  change  in  dse,  which  causes  an  increase  in  gav,  than  it  is  for  positive.  Hence  when 

2  choosing  the  bias  potential  F*  of  HARP  for  a  desired  £&>,  the  smallest  dse  for  the  entire  detector  area 

3  should  be  used,  so  that  the  variation  in  thickness  is  positive  from  the  reference  thickness.  If  we 

4  assume  the  smallest  dse  =  8  pm,  the  percentage  changes  in  gav  due  to  an  increase  in  dse  of  1  %  are  19 

5  %,  23  %  and  26  %,  respectively,  for  the  three  F*  settings  of  840  V,  864  V  and  880  V.  The  change  in 

6  gaV  increases  to  33  %,  40  %,  and  44  %,  respectively,  when  the  variation  in  dse  increases  to  2  %.  This 

7  range  of  gav  non-uniformity  is  certainly  within  the  capability  of  gain  correction  and  the  dynamic 

8  range  of  the  electronic  circuits  of  FPI.  With  further  increase  in  dse  to  4  %,  the  value  for  gav  will 

9  decrease  from  the  original  setting  of  46,  25  and  12  to  15,  10  and  6,  respectively,  which  is  a  gain  non- 

10  uniformity  of  2  to  3  folds.  Although  the  variation  in  gav  is  significant,  it  is  not  expected  to  cause 

1 1  significant  degradation  in  the  DQE  performance,  as  can  be  seen  in  the  example  of  gav  =  46  and  12 

12  shown  in  Figure  7  (a).  However  this  will  pose  a  greater  challenge  for  the  uniformity  correction 

13  algorithm.  Hence  we  can  conclude  that  a  thickness  non-uniformity  of  2  %  is  desirable  for  a  HARP 

14  layer,  and  4  %  non-uniformity  could  still  possibly  be  tolerated. 

15 

16  V.  CONCLUSIONS 

17  The  feasibility  of  a  new  concept  of  indirect  conversion  flat-panel  detector  with  avalanche  gain  is 

18  being  investigated.  The  avalanche  gain  and  quantum  efficiency  of  an  8  pm  thick  HARP  layer  were 

19  investigated  experimentally  as  function  of  Ese •  These  characteristics  of  HARP  were  applied  to  a 

20  cascaded  linear  system  to  determine  the  optimal  operating  conditions  required  for  x-ray  imaging. 

21  Our  results  showed  that  by  varying  the  avalanche  gain  between  2  and  46,  x-ray  quantum  noise 

22  limited  performance  and  a  wide  dynamic  range  can  be  obtained  for  both  R/F  and  mammographic 

23  imaging  applications.  Direct  x-ray  interaction  in  HARP  was  shown  not  to  affect  the  imaging 


21 


1  performance.  Due  to  the  rapid  increase  in  gav  as  a  function  of  Ese,  a  small  variation  in  dse  can  cause 

2  significant  non-uniformity  in  avalanche  gain.  It  is  necessary  to  keep  the  thickness  uniformity  of  a 

3  HARP  layer  within  4  %,  which  could  maintain  a  high  DQE  across  the  entire  detector  and  make  the 

4  variation  in  gav  correctable  using  gain  calibration. 
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Table  I:  SHARP -AMFPI  detector  design  parameters  and  operating  conditions  chosen  for  different  x 
ray  imaging  applications. 


Mammography  detector 

R/F  detector 

tomosynthesis 

screening 

Fluoroscopy 

Pixel  pitch  av  (pm) 

70 

200 

X-ray  spectrum 

28  kVp  Mo/Mo,  4cm  PMMA 

RQA5 

Electronic  noise 

2000  electrons 

1500  electrons 

Csl  (75%  packing) 

150  pm  HR 

600  pm  HL 

W  for  CsI(eV) 

30 

18 

a 

0.8 

0.8 

Minimum  exposure 

0.1  mR 

1  mR 

0.1  pR 

30  pR 

Ese  (V/pm) 

110 

105 

110 

95 

0 

0.36 

0.35 

0.36 

0.31 

Rav 

46 

12 

46 

2.3 

Figure  Caption 

Figure  1 :  (a)  Diagram  showing  the  concept  of  SHARP- AMFPI,  where  an  avalanche  photoconductor  is 
used  to  detector  light  photons  generated  by  Csl.  (b)  Side  view  showing  the  multilayer  structure  of  the 
detector  (thickness  not  drawn  to  scale). 

Figure  2:  Experimental  apparatus  for  measuring  the  photosensitivity  of  HARP. 

Figure  3:  Measured  effective  quantum  efficiency  q*  as  a  function  of  Ese  for  an  8  pm  thick  HARP  layer. 
Figure  4:  The  optical  quantum  efficiency  of  a-Se  calculated  using  Onsager  theory  and  r0  =  1.7  nm, 
compared  with  the  measured  effective  quantum  efficiency  q  at  X  =  540  nm. 

Figure  5:  (a)  Solid  circles  are  the  gav  for  the  8  pm  HARP  layer  calculated  by  dividing  the  measured  q* 
at  540  nm  by  the  optical  q  predicted  by  the  Onsager  theory  using  r0  =  1.7  nm.  The  solid  line  shows  the 
gav  calculated  using  the  fitted  Pi  and  p2  values,  (b)  The  plot  of  ln(ln(gflV))  as  a  function  of  MEse,  and  the 
best  linear  fit  for  the  data. 

Figure  6:  Flow  diagram  showing  the  serial  and  parallel  stages  of  the  cascaded  linear  system  model  for 
SHARP-AMFPI. 

Figure  7:  Calculated  DQE(/)  using  the  detector  parameters  and  operating  conditions  shown  in  Table  I 
for  a  fluoroscopy  detector  at  the  x-ray  exposure  level  of:  (a)  0. 1  pR;  and  (b)  30  pR. 

Figure  8:  The  calculated  image  charge  on  each  pixel  of  the  detector  and  the  corresponding  avalanche 
gain  as  a  function  of  x-ray  exposure  for  the  detector  operational  setting  of:  (a)  fluoroscopy;  and  (b) 
radiography. 

Figure  9:  Calculated  DQE(/)  using  the  detector  parameters  and  operating  conditions  shown  in  Table  I 
for  mammography:  (a)  tomosynthesis  with  minimum  exposure  of  0.1  mR;  and  (b)  screening  with 
minimum  exposure  of  1  mR. 

Figure  10:  The  calculated  image  charge  on  each  pixel  of  the  mammography  detector  and  the 
corresponding  avalanche  gain  as  a  function  of  x-ray  exposure  with  or  without  the  addition  of  a  pixel 


storage  capacitor.  The  results  are  for  detector  operating  conditions  chosen  for:  (a)  tomosynthesis  image 
acquisition;  and  (b)  screening  mammography. 

Figure  1 1 :  Diagram  showing  the  depth  dependent  gain  for  direct  x-ray  interaction  in  HARP.  For  the  x- 
rays  absorbed  in  Csl,  the  light  photons  created  enter  the  top  interface  of  HARP,  where  holes  are 
generated  and  undergo  avalanche  multiplication  for  the  entire  thickness  of  HARP  dse-  For  the  x-rays 
interacting  directly  with  HARP  at  depth  x  from  the  top  interface,  the  holes  generated  undergo 
avalanche  for  the  thickness  of  ( dse  -  x),  which  causes  depth  dependence  avalanche  gain  variation. 

Figure  12:  Comparison  of  the  presampling  signal  and  NPS  (before  pixel  aperture  function)  due  to  x-ray 
absorbed  in  Csl  and  direct  x-ray  interaction  in  HARP:  (a)  signal  spectra,  where  the  signal  due  to  direct 
interaction  is  negligible;  and  (b)  NPS,  where  the  NPS  due  to  direct  interaction  under  the  Nyquist 
frequency  of  2.5  cycles/mm  is  negligible. 

Figure  13:  Comparison  between  the  DQE  of  SHARP-AMFPI  with  and  without  consideration  of  direct 
x-ray  interaction  in  HARP. 

Figure  14:  The  avalanche  gaingav  calculated  as  a  function  of  the  thickness  of  HARP,  which  varies 
around  8  pm,  under  a  constant  bias  potential  V of  840,  864  and  880  V. 
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